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SECTION  1 


INTRODUCTION 

Continuous  flow  combustion  systems  are  both  turbulent  and  strongly 
backmixed.  Such  flows  are  classified  as  complex  because  of  the  complicated 
fluid  dynamics  in  contrast  to  simple  flows  which,  while  turbulent  as  well  as 
laminar,  do  not  have  zones  of  recirculation  (Launder,  1978).  Examples  of 
practical  continuous  flow  combustion  systems  of  direct  interest  to  the  Air 
Force  are  gas  turbine  and  ramjet  combustors. 

An  increased  understanding  of  complex  flows  is  required  as  gas  turbine 
and  ramjet  combustors  are  pushed  to  higher  levels  of  performance,  and  as  the 
fuel  available  for  gas  turbine  engines  decreases  in  quality.  This  need, 
somewhat  fortuitously,  coincides  with  the  advances  in  laser  diagnostics  and 
modeling  that  will  be  necessary  to  provide  the  needed  insight. 

Experiments  are  required  in  laboratory  systems  that  provide  spatially 
resolved  measurements  of  velocity,  temperature,  and  concentration  in  complex 
flows  representative  of  both  premixed  and  non-premixed  systems.  Mathematical 
modeling  is  required  in  order  to  plan  and  Interpret  experimental  research  and, 
ultimately,  to  design  hardware. 

The  goal  of  the  present  AFOSR  grant  at  the  UCI  Combustion  Laboratory  was 
to  develop  laboratory  model  combustors  and  experimental  methodology  suitable 
for  the  acquisition  of  the  desired  Information.  The  objectives  of  the  program 
were  as  follows: 

(1)  To  establish  and  evaluate  laboratory  model  combustors  for  both 
premixed  and  non-premixed  fuel/air  injection  suitable  for  studies  of 
complex  flows, 

(2)  To  acquire  spatially-resolved  flowfleld  data,  using  state-of-the-art 
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diagnostics  for  the  purposes  of  model  verification  and  providing 
Insight  Into  the  transport  of  mass,  momentum,  and  energy  In  complex 
flows, 

(3)  To  conduct  supplementary  studies  to  support  the  development  of 
methodologies  required  for  measurements  in  complex  flows. 

To  meet  these  objectives,  the  program  was  divided  in  a  task  organization 
combusting  of  three  elements: 

•  Element  A:  Model  Laboratory  Combustor  Development  and  Evaluation 

•  Element  B:  Physical  and  Chemical  Processes  in  Complex  Flows 

•  Element  C;  Supplemental  Studies 

The  experimental  facilities  developed  to  support  the  research  are 
described  in  the  following  section.  In  Section  3.0,  the  modeling  procedure 
used  in  the  present  program  is  described.  The  results  are  presented  in 
Section  4.0  for  each  of  the  elements  delineated  above.  Finally,  a  summary  is 
provided  in  Section  5.0. 
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SECTION  2 

EXPERIMENTAL  FACILITIES 
2.1  Laboratory  Combustor  Facility 

The  Laboratory  Combustor  Facility  (LCF)  provides  air  flow  sufficient  to 
produce  a  reference  velocity  of  15.24  mps  (50  ft/sec)  In  a  80  ms  (3-lnch) 
duct,  and  flow  rates  of  gaseous  fuels  to  yield  overall  equivalence  ratios 
exceeding  2.0.  A  separate  fuel  and  air  supply  system,  Independent  of  the  main 
flow.  Is  available  for  operating  the  jet  flow  In  the  opposed-jet  combustor 
module.  Details  of  the  test  stand  are  available  (Peterson  and  Himes,  1978). 

The  facility  Is  designed  to  accommodate  a  variety  of  combustor 
configurations.  In  the  present  research  program,  the  following  four  combustor 
modules  were  evaluated: 

Combustor  Module 
Opposed-Jet 
Centerbody 

Centerbody  with  Swirl 
Dilute  Swirl 

A  schematic  of  each  module  Is  presented  In  Figure  1. 

Opposed  Jet  Combustor  Module.  The  opposed  jet  combustor  (OJC)  module 
consists  of  a  51  nm  (2-lnch)  Inside  diameter  cylindrical  Vycor  duct  containing 
an  opposing  axlsymmetrlc  jet  (Figure  la).  The  Inner  diameter  of  the  tube  from 
which  the  jet  Issues  Is  1.3  mm  (0.052-lnch).  The  outer  diameter  of  the  tube 
Is  6.4  mn  (0.25-lnch).  The  space  between  the  Inner  and  outer  walls  Is  hollow 
except  for  one  partitioning  barrier  to  provide  for  two-pass  water  cooling  (In 
the  case  of  HOT,  reacting  flow  studies)  or  water  heating  (In  the  case  of 


Fuel/ Air  Injection 
Premixed 
Non-Premlxed 
Non-Premlxed 
Non-Premlxed 


•EXIT  PLANE* 


Dilute  Swirl  Combustor 


32cm 


Figure  1.  (concluded) 


HEATED  jet,  non-reacting  flow  studies). 

The  Vycor  duct  Is  butted  against  and  attached  to  an  equal  diameter 
stainless  steel  duct  18  duct  diameters  (l.e.,  92  cm)  In  length.  The  purpose 
of  the  stainless  steel  duct  Is  to  provide  some  relaxation  of  the  flow  before 
entering  the  Vycor  duct.  The  Vycor  duct  Is  nominally  9  duct  diameters  In 
length  (46  cm).  For  purposes  of  measurement,  the  Inlet  plane  (l.e.,  x/d  ■  0) 
Is  arbitrarily  assigned  an  axial  location  4.25  duct  diameters  downstream  of 
the  joining  point  of  the  Vycor  and  stainless  steel  duct.  This  places  the  exit 
plane  of  the  duct  at  x/d  >  4.71  (l.e.,  240  nm  from  the  Inlet  plane).  The 
Inlet  plane  of  the  opposing  jet  Is  arbitrarily  selected  to  be  x/d  -  3.54 
(180  ran  from  the  Inlet  plane). 

Details  of  the  experimental  setup  are  described  by  Peterson  and  Himes 
(1978),  and  Wuerer  (1978). 

Centerbody  (Non-Swlrl)  Module.  The  centerbody  combustor  (CBC)  module  is 
a  l/5th  scaled  version  of  the  combustor  developed  and  used  at  the  Air  Force 
Wright  Aeronautical  Laboratory  (Roquemore,  et  al.,  1980).  The  CBC  consists  of 
a  51  mm  (2-lnch)  ID  cylindrical  pyrex  tube  (3  ran  wall)  that  houses  an  Internal 
centerbody  of  30.5  ran  (1.2-inch)  OD  around  which  the  annular  air  can  be 
swirled  (Figure  lb).  Fuel  (Propane  for  reacting  HOT  cases;  CO2  for  non¬ 
reacting  COLO  and  HEATED  cases)  Is  Introduced  through  a  jet  of  1.3  mm  (0.052- 
Inch)  diameter  centrally  located  In  the  centerbody.  The  test  section  Is 
arbitrarily  selected  to  be  at  20  cm  (7.87-inch)  as  measured  from  the 
centerbody  face.  The  overall  length  of  the  tu:.-  Is  70  cm  ( 27. 55-lnch). 

Details  of  the  experimental  setup  are  available  (Ikloka,  1982). 

Centerbody  (Swirl)  Module.  Swirl  was  added  to  the  CBC  In  the  present 
study  to  evaluate  the  performance  of  the  module  with  a  swirl-induced 
recirculation  zone  (Figure  Ic).  The  swirl  was  Induced  In  the  annular  air  by 
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placement  of  45  degree  straight  swirl  vanes  at  a  location  57  mm  (2.24-lnch) 
upstream  of  the  centerbody  face. 

Details  of  the  experimental  setup  are  available  In  Ikloka  (1982). 

Dilute  Swirl  Module.  The  Dilute  Swirl  Combustor  (DSC)  module  consists  of 
a  80  mm  (3-lnch)  ID  duct  (Figure  Id)  within  which  is  located  a  small  central 
body  just  sufficient  to  house  a  small  fuel  nozzle.  Immediately  surrounding 
the  central  body  Is  a  set  of  swirl  vanes.  Surrounding  the  swirl  vanes  is  non- 
swlrled  dilution  air.  Dilution  and  swirl  air  are  metered  separately.  The 
dilution  (l.e.,  non-swlrllng)  air  Is  introduced  through  a  flow  stralghtener  in 
the  outer  annulus.  The  swirl  air  passes  through  a  set  of  swirl  vanes  which 
impart  a  solid  body  angle  of  turn  to  the  flow,  60°  In  the  present  case. 

Either  liquid  or  gaseous  fuels  may  be  Introduced  through  a  nozzle  at  the  end 
of  the  central  tube.  The  exit  plane  of  the  fuel  injector  Is  set  at  the  same 
axial  location  as  the  exit  plane  of  the  swirl  air  to  provide  a  clean,  well 
defined  boundary  condition  for  modeling.  The  overall  combustor  tube  length 
may  be  varied. 

Details  of  the  experimental  setup  are  available  (Brum,  1983). 

2.2  Hydrodynamic  Facility 

The  Rydrodjmamlc  Facility  (HDF)  Is  a  continuous  flow,  closed-loop  system 
with  water  circulated  by  a  centrifugal  pump.  The  test  section  consists  of  two 
76  mi  (3-lnch)  ID  tubes,  380  cm  (12.5-foot)  and  61  cm  (24-lnch)  In  length 
respectively,  mounted  horizontally.  Combustor  modules,  such  as  the  opposed 
jet  and  centerbody  combustor  modules,  are  Installed  between  the  two  lengths  of 
pipe.  In  the  case  of  the  OJC,  the  380  cm  long  tube  serves  as  the  test 
section.  In  the  case  of  the  CBC,  the  61  cm  long  tube  serves  as  the  test 
section. 

The  opposed  jet  configuration  adopted  for  the  present  study  (Figure  2a) 
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employs  an  0.198  cm  (0.078-inch)  inner  diameter  (ID)  jet  Cube,  having  an  0.952 
cm  (0.375-inch)  outer  diameter  (OD),  and  made  of  type  316  stainless  steel. 

The  jet  tube  is  positioned  concentrically  in  a  76  im  (3-inch)  ID  Pyrex  glass 
duct  which  comprises  Che  test  section  of  the  RDF.  The  Jet  Cube  extends  23.5 
cm  (9.25-inch)  into  the  test  field.  The  opposed  jet  has  a  mass  flow  rate  of 
approximately  1  i/m  and  opposes  mainstream  flows  of  87  i/m  and  174  i/m  in  two 
different  test  cases. 

The  centerbody  configuration  (Figure  2b)  consists  of  a  4.41  cm  (1.74- 
inch)  OD  cylindrical  PVC  bluff  body,  12.7  cm  (5-lnch)  in  length,  with  a  flat 
end  face.  The  centerbody  is  situated  concentrically  in  Che  76  an  (3-inch)  ID 
glass  duct  of  the  HDF.  An  axisymmetric  jet  issues  from  an  0.198  cm  (0.078- 
inch)  diameter  opening  in  the  center  of  Che  bluff  body  face  and  has  a  flow 
rate  equal  to  0.64  percent  of  the  total  flow  (46.3  gpm)  of  Che  system. 

Details  of  the  experimental  setup  are  available  in  Arcoumanis  (1980)  and 
Grudovich  (1983). 

2.3  Laser  Anemometer  (LA) 

Two  laser  anemometers  were  developed  and  employed  in  the  present 
program:  a  single  component  and  a  two  component.  The  former  is  used 
exclusively  on  the  HDF  whereas  the  latter  supports  the  LCF. 

Single  Component  LA.  The  single  component  LA,  shown  in  Figure  3a,  is  a 
dual  beam  anemometer.  The  laser  beam  is  split  into  two  equal  intensity  parts 
which  are  made  to  intersect  by  means  of  a  system  of  plano-convex  lenses.  As  a 
laser  source,  a  5mW  Hellum-Heon  laser  (Spectra-Physics  model  #120)  is  used 
which  operates  in  the  TEM  mode,  Che  lowest  order  transverse  mode.  The  0.8  mm 
beam  diameter  is  focused  by  means  of  lens  L^  ^^1  ~  on),  which  is  mounted 
on  the  laser  head,  onto  a  radial  diffraction  gracing.  The  Technlsch  Physische 
Dienst  (TPD)  rotating  grating  functions  both  as  a  beam  splitter  and  as  a 
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a)  Single-Component 


b)  Two-Component 


MAINSTREAM 


Figure  3.  Laser  Anemometer  (LA) 


frequency  shifter.  As  a  beam  splitter,  it  diffracts  the  Incident  laser  beam 
Into  two  equal-intensity  beams  Independent  of  the  state  of  polarization  of  the 
incident  light.  On  the  other  hand,  as  a  frequency  shifter,  it  provides  a 
frequency  shift  In  the  two  diffracted  beams,  thus  causing  a  moving  fringe 
pattern  In  a  direction  which  depends  on  the  direction  of  the  disc  rotation. 

The  laser  anemometer  setup  Is  mounted  on  an  optical  bench  and  the  whole 
system  Is  traversed  In  the  horizontal  plane  along  the  axis  of  the  test  section 
and  perpendicular  to  this  axis.  The  optical  arrangement  results  In  a  fringe 
volume  diameter  of  0.76  mn  and  length  of  0.58  mm.  The  scattered  light  from 
the  particles  which  are  traversing  the  fringes  Is  detected  by  a  DISA 
photomultiplier  (Model  <!>55L10). 

The  signal  is  processed  by  either  a  DISA  Model  55L20  Frequency  Tracker  or 
a  SDL  Model  VP  1001  Counter  Processor.  Details  of  the  experimental  setup  are 
provided  in  Acroumanls  (1980)  and  Grudovich  (1983). 

Two-Component  LA.  The  two-color  laser  anemometry  (LA)  system  is  shown  In 
Figure  3b.  The  beam  from  a  200  mW  Argon-ion  laser  (Lexel  Model  75)  is 
collimated  and  passed  through  a  prism  to  separate  the  various  wavelengths. 

The  two  most  Intense  beams,  green  (514  nm)  and  blue  (488  nm)  are  each  passed 
through  a  series  of  optics  In  which  they  are  polarized  and  split  Into  two 
beams  of  equal  Intensity  50  mm  apart.  An  upstream  40  MHz  frequency  shift  (TSI 
model  915  Bragg  Cell)  Is  applied  to  one  of  each  pair  of  beams  In  order  to 
avoid  directional  ambiguity  that  would  otherwise  result  from  the  highly 
turbulent  recirculating  flow.  The  four  beams  (blue  pair  In  the  vertical  plane 
and  green  pair  In  the  horizontal  plane)  are  then  focused  through  a  250  mm  lens 
to  a  common  point  within  the  test  section.  This  results  In  set  of 
perpendicular  Interference  fringes  spaced  at  2.6  ms  the  green  beams 
(vertical  fringes)  and  2.5  izn  for  the  blue  beams  (horizontal  fringes)  which 
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are  responsive  to  the  axial  and  tangential  velocity  components  respectively. 

Receiving  optics  consist  of  a  120  nn  lens  focused  onto  a  0.25  mm  diameter 
photomultiplier  tube  aperture  (via  an  appropriate  dlchromate  filter  to 
selectively  pass  either  the  blue  or  green  light).  These  optics  are  placed  at 
an  angle  of  20**  off  direct  forward  scatter  which  results  In  a  probe  volume  of 

0.022  mm  and  cross-sectional  area  perpendicular  to  the  axis  of  measurement  of 

2 

0.10  mm  .  However,  due  to  the  requirement  Imposed  by  the  processing 
electronics  that  both  axial  (u)  and  tangential  (w)  velocity  components  be 
obtained  simultaneously,  the  effective  probe  cross-section  Is  much  less 
(approximately  0.03  mm  ).  The  transmitting  and  receiving  optics  are  mounted 
on  an  optical  bench  capable  of  placing  the  measurement  volume  at  points 
throughout  the  stationary  combustor  test  section. 

A  special  electronic  Interface  was  built  to  Interface  the  output  of  the 
two  counter  process  channels  (u,w)  directly  to  a  DEC  PDF  11/23  computer 
system.  This  Interface  Identifies  whether  or  not  the  u  and  w  events  occur 
within  a  certain  aperature  time  of  each  other  (normally  ^  50  us)*  H  so, 
they  are  considered  simultaneous,  stored,  and  then  multiplexed  Into  the  11/23 
via  a  parallel  Interface.  Once  the  Interface  verifies  that  the  11/23  has  read 
the  data.  It  simultaneously  resets  both  processor  channels.  The  key  feature 
of  this  system  Is  that  It  permits  a  direct  and  Instantaneous  measurement  of 
correlations  such  as  u'w’ .  An  aperature  time  of  50  us  was  selected  since,  at 
the  maximum  bulk  velocity  measured  (15  m/s),  an  equivalent  spatial  resolution 
of  less  than  0.75  mm  Is  obtained. 

The  11/23  computer  Is  equipped  with  an  Internal  clock  having  a  resolution 
of  100  us  which  Is  Initiated  at  the  beginning  of  each  run  cycle.  As  the  u,  w 
data  are  received,  the  time  of  event  t  Is  combined  with  the  raw  data  (u,w,t), 
and  the  data  are  permanently  stored  in  an  archival  fashion  for  future 
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reference  and  analysis. 

Details  of  the  system  are  provided  In  Brum  (1983). 

Seeding.  Natural  seeding  Is  used  In  the  HDF.  In  the  LCF,  the  main  and 
fuel  jet  flows  are  seeded  Independently,  but  to  the  same  levels  of 
concentration,  with  1  ym  alumina  particles.  A  liquid  suspension  atomization 
seeding  technique  Is  employed.  Signal  validation  was  obtained  using  two 
counter  processors  (Macrodyne  Model  2098).  This  seeding  technique,  developed 
under  the  present  program.  Is  described  In  Section  4.3.1. 

2.4  Temperature  Measurements 

Mean  temperatures  were  obtained  with  a  robust  thermocouple  probe  and, 
under  limited  combustor  operating  conditions,  fluctuating  and  mean 
temperatures  were  obtained  with  a  fine-wire  thermocouple  probe. 

Robust  Thermocouple  Probe.  The  robust  probe  consists  of  an  unshielded, 
fine  wire,  platinum/ platinum  13%  rhodium  thermocouple  mounted  on  a  micrometer 
traverse  jig.  A  Doric  digital  pyrometer  (Model  DS-500)  is  used  to  record  the 
thermocouple  output. 

Details  of  the  probe  are  provided  in  McDannel  (1979). 

Fine-Wire  Thermocouple  Probe.  A  fine-wire  probe  was  built  under  a 
complementary  National  Science  Foundation  (NSF)  Grant  (CPE-8013742)  and  used 
In  conjunction  with  the  two-component  LA  to  acquire  measurements  of  axial  and 
azimuthal  heat  flux.  Platinum/platinum  10%  rhodium  thermocouple  wire  of  three 
diameters  (25  gm,  SO  ym,  and  125  ym)  are  used  In  this  work.  The  thermocouple 
junction  Is  formed  by  overlapping  and  spotwelding  the  two  small  wires  which 
yields  a  junction  with  a  characteristic  length  less  than  40  ym  when 
25  m  wire  Is  used.  The  small  diameter  wire  Is  then  gas  welded  to  larger 
support  wires  (of  250  in  and  500  ym  diameter)  of  the  same  material.  The' 
support  wires  are  cemented  In  a  0.159  cm  (0.063-lnch)  or  0.318  cm  (0.125-lnch) 
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OD  alumina  tube  which  is  in  turn  placed  in  Inconel  tubing  of  various  diameters 
and  shapes.  The  various  probe  configurations  are  tested  and  the  one  that 
minimizes  perturbation  effects  at  a  specific  measurement  location  is  used. 
Figure  4a  shows  the  dimensions  of  the  probe  configurations  used  during  these 
experiments.  The  probe  is  held  in  a  mechanical  traverse  which  moves  with  the 
LA  optical  table  and  is  inserted  into  the  flow  through  the  exhaust  plane  of 
the  combustor. 

The  thermocouple  frequency  response  necessary  to  make  accurate  time- 
resolved  measurements  can  be  estimated  from  the  physical  dimensions  of  the 
combustor  and  the  bulk  mean  velocity.  If  the  length  scale  of  the  energy- 
containing  structure  is  taken  to  be  one-half  the  combustor  radius,  the 
corresponding  frequency  for  a  mean  velocity  of  15  m/s  is  about  750  Hz. 
Therefore,  in  order  to  accurately  measure  the  root-mean-square  temperature, 
the  sensor  frequency  response  must  extend  to  at  least  two  or  three  times  the 
frequency  of  the  energy-containing  structures  or  to  approximately  2  kHz,  The 
uncompensated  frequency  response  of  the  25  pm  sensor  at  this  mean  velocity  is 
about  15  Hz.  Thus,  the  frequency  response  must  be  Increased  by  a  factor  of 
130  and  is  accomplished  using  electronic  compensation. 

The  compensation  method  used  in  this  study  is  similar  to  the  one  used  by 
Lockwood  and  Monelb  (1981).  A  block  diagram  of  the  thermocouple  signal 
processing  electronics,  shown  in  Figure  4b,  can  be  divided  into  the  following 
subcircuits:  1)  a  relay,  a  source  of  heating  current,  and  a  square  wave 
generator  to  control  the  relay,  2)  a  low  noise,  low  drift,  differential  input 
preamplified  (AD  521;  gain  of  100),  3)  a  compensation  circuit,  and  4)  a  gain 
and  bias  amplifier  to  match  the  dynamic  range  of  the  temperature  signal  with 
that  of  the  analog  to  digital  converter. 

The  compensation  subcircuit  consists  of  a  differentiator,  an  amplifier. 


14 


a)  Probe  Configurations 
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b)  Electronics 
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Figure  4.  Fine-Wire  Thermocouple 


and  a  sunualng  amplifier.  The  operation  [1  +  t  (d/dt)]E,  where  x  is  the  mean 
time  constant,  is  performed  on  the  thermocouple  voltage,  E,  to  increase  the 
inherently  low  (~  15  Hz)  frequency  response  of  the  sensor. 

Details  of  the  experimental  setup  are  available  (Seiler,  1983a). 

2.5  Emissions  Measurement 

Combustion  gases  are  sampled  continuously  using  a  hot  (60“C/140®F) 
watercooled,  6.35  mm  (0.25-lnch)  OD  by  1.32  mm  (0.052-inch)  ID,  stainless 
probe  with  a  constant  area  inlet  bounded  by  a  cooled,  streamlined  tip.  The 
sample  probe  is  mounted  on  a  micrometer  traverse  Jig  for  sampling.  The  gas 
sample  flows  to  the  analysis  instruments  through  heated  (65®C/ 150®F) ,  6.35  mm 
(0.25-inch)  OD  Teflon  tubing. 

Analysis  of  the  sample  gas  is  performed  using  a  packaged  emission 
analysis  system  (Scott  Research,  Model  113).  Passage  of  the  gas  through  an 
ice  bath  allows  all  concentration  measurements  to  be  taken  on  a  dry  basis. 
Nitrogen  oxides  (NO,  NO^^)  concentrations  are  measured  using  a 
chemiluminescence  analyzer  (Scott  Research,  Model  125).  Nondispersive 
Infrared  analyzers  (Beckman,  Models  315B  and  315BL  respectively)  are  used  to 
measure  carbon  dioxide  (CO2)  and  carbon  monoxide  (CO)  concentrations.  Total 
hydrocarbon  (HC)  concentration  is  measured  by  a  flame  ionization  detector 
(Scott  Research,  Model  215).  Oxygen  (O2)  concentration  is  measured  by  a 
paramagnetic  analyzer  (Scott  Research,  Model  150). 

Details  of  the  experiment  setup  are  provided  in  McDannel  (1979). 


SECTION  3 


MODELING 

3.1  General  Methodology 

The  computational  procedure  used  to  obtain  the  simultaneous  solution  of 
the  governing  elliptic  partial  differential  equations  Is  the  TEACH  method 
(Gosman  and  Pun,  1974)  which  has  been  expanded  (1)  for  the  opposed  jet, 
centerbody,  and  dilute  swirl  configurations,  and  (2)  to  Include  species 
conservation,  swirl,  and  reaction. 

The  governing  equations  are  cast  Into  a  generalized  equation  of  the  form: 


h  7  If  ^^*,eff 


r  3r  ^*^^4.,eff  Ir^  ^(| 


(1) 


where : 

(p  •  dependent  variable 

r  cf  "  effective  exchange  coefficient 

(>,err 

S.  *  source  and/or  sink  terms 

The  boundary  conditions  are  specified  along  the  four  boundaries  of  the 
solution  domain.  Where  applicable,  experimental  data  are  used  to  specify  the 
boundary  conditions. 

3.2  Turbulence  Models 

The  form  of  the  governing  equations  Is  based  upon  the  gradient  transport 
hypothesis.  Closure  of  the  equations  requires  the  development  of  time  mean  of 
"effective”  transport  coefficients.  The  effective  viscosity  defined 


as  follows: 


(2) 


‘^eff 


where: 

■  absolute  viscosity  of  fluid 

■  eddy  (or  apparent)  viscosity  due  to  turbulent  exchange 

The  turbulent  transport  of  scalar  quantities  is  related  to  the  eddy  viscosity 
through  the  appropriate  Schmidt  numbers  as  follows. 


^♦.eff 

where : 


<j  <J 

♦.t 


laminar  Schmidt  Number 
turbulent  Schmidt  Number 


(3) 


Two  models  were  employed  for  the  turbulent  viscosity  (y^)  ,  an  algebraic 
model  and  a  two-equation  k-e  model. 

Algebraic  Hodel.  The  employment  of  an  algebraic  model  simplifies  the 
closure  of  the  governing  equations  be  eliminating  the  need  for  the  equations 
for  k  and  e.  The  use  of  algebraic  models  for  recirculating  flowfields  has 
shown  little  evidence  of  success.  However,  their  application  is  useful  in 
revealing  solution  sensitivities  to  viscosity  distribution.  The  algebraic 
model  utilized  in  the  present  Investigation  was  one  proposed  for  the 
preliminary  demonstration  of  numerical  combustion  models. 
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where : 


K  *  empirical  coefficient 
D  >  diameter  of  combustor 
L  >  length  of  combustor 
p  >  local  density 
mU^  ■  inlet  kinetic  energy 

As  can  be  seen,  the  variation  in  is  dependent  only  upon  local  density  for  a 
given  geometry  and  inlet  flow  rate.  For  the  Isothermal  flow  of  a  mixture  of 
two  gases  with  assumed  Identical  density  and  laminar  transport  properties,  the 
case  of  the  algebraic  viscosity  model  becomes  extremely  simple  (l.e.,  the 
effective  viscosity,  Ugf£i  la  constant  through  the  flowfleld). 

Two-Equation  k-e  Model .  The  two-equation  k-e  model  is  the  most  commonly 
employed  viscosity  model  for  combustor  flow  analyses.  The  two-equation  model 
accounts  for  the  transport  effects  of  two  turbulence  characteristics  which  may 
be  taken  as  typical  time  scales  and  length  scales  of  the  energy  containing 
motions.  This  is  accomplished  by  calculating  the  local  levels  of  turbulence 
kinetic  energy  (k)  and  the  kinetic  energy  dissipation  rate  (e).  The 
representative  length  and  time  scales  are: 


I 


(5) 


and 


e 


(6) 


The  eddy  viscosity  is  then  calculated  from  the  algebraic  Prandtl-Kolmogorov 


formulation,  l.e 


(7) 


where  Is  an  empirical  coefficient. 

Detailed  information  on  the  modeling  in  general  is  provided  elsewhere 
(Wuerer,  1978;  Brum,  1983,  and  Grudovlch,  1983). 
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SECTION  4 


RESULTS 

A  summary  of  the  results  are  provided  In  this  section  for  each  of  the 
Elements,  A,  B,  and  C.  Details  are  provided  In  papers  presented, 
publications,  and  both  Master  of  Science  Theses  and  Doctor  of  Philosophy 
Dissertations  completed  under  the  present  grant. 

4.1  Element  A;  Model  Laboratory  Combustor  Developments  and  Evaluation 

The  Opposed  Jet  Combustor  (OJC)  module  was  developed  under  a  previous 
AFOSR  grant  (74-2710)  as  a  candidate  configuration  for  a  premixed  model 
complex  flow  laboratory  combustor.  In  the  present  program,  attention  was 
directed  to  candidate  configurations  for  a  non-premixed  model  complex  flow 
laboratory  combustor  using  the  following  criteria  for  the  purposes  of 
evaluation: 

•  allow  good  optical  access  for  advanced  laser  diagnostics, 

e  be  stabilized  by  a  zone  of  strong  backmixlng, 

•  be  capable  of  swirl  Inlet  conditions, 

•  be  axlsymmetrlc  to  simplify  modeling, 

•  have  clean,  well  defined  boundary  conditions, 

•  be  capable  of  reacting  and  non-reacting,  non-premlxed  (multiple 
Inlet  stream)  operation, 

•  allow  sampling  probe. 

As  a  point  of  departure,  a  1/5  scale  model  (Figure  lb)  of  a  geometry 
Introduced  for  diagnostics  and  modeling  development  at  the  Air  Force  Wright 
Aeronautical  Laboratory  (Roquemore,  et  al.,  1980)  was  evaluated.  The  scaled 
version  of  the  Centerbody  Combustor  (CBC)  consists  of  a  51  mm  ID  cylindrical 
pyrex  tube  (3  nm  wall)  that  houses  an  Internal  centerbody  of  30.5  am  O.D. 


around  which  the  air  is  Injected  in  an  annular,  non-swirling  jet. 

The  CBC  was  evaluated  in  the  present  study,  not  only  for  the  nominal 
operating  mode  without  swirl,  but  with  a  swirl-induced  annular  jet  as  well 
(Figure  Ic).  Fuel  (gaseous  propane  for  reacting  cases,  CO2  for  non-reacting 
cases)  was  Introduced  through  the  jet  centrally  located  in  the  centerbody. 

In  addition  to  the  CBC  geometry,  an  independent  ’’Dilute  Swirl  Combustor" 
(DSC)  geometry  (Figure  Id)  was  designed  to  satisfy  the  established  criteria 
and  address  limitations  observed  with  the  no-swirl  and  swirl  operation  of  the 
CBC.  In  addition  to  propane  for  the  reacting  cases  and  CO2  for  the  non¬ 
reacting  cases,  the  DSC  was  also  operated  with  liquid  Isooctane  injected 
through  either  a  twin-fluid  or  pressure  atomizer. 

Both  the  CBC  and  DSC  modules  were  operated  at  atmospheric  pressure,  with 
ambient  inlet  temperature.  Axial  mean  and  rms  velocity  measurements  were  made 
using  the  two  component  laser  anemometry  system  (Figure  3b). 

High-speed  photography  was  used  to  qualitatively  document  the  flame 
structure  and  investigate  dynamic  aspects  of  the  flame.  The  camera  (Fastax 
Model  WF3)  was  run  at  a  frame  speed  of  approximately  1000  frames/second. 

Black  and  white  film  (Kodak  7277,  16  mm,  400  ASA)  was  used  and  pushed  2  stops 
during  processing. 

4.1.1  No-Swirl  CBC 

Swirl  is  one  of  the  primary  features  required  by  a  test  bed  geometry  to 
be  useful  as  a  supporting  data  case  for  the  identified  modeling  efforts.  The 
no-swirl  CBC  was  investigated  as  a  point  of  reference  and  because  it  is  the 
configuration  used  at  the  Air  Force  Wright  Aeronautical  Laboratory  (AFWAL). 

Laser  anemometry  measurements  were  made  to: 

(1)  Identify  the  basic  flow  patterns  including  the  extent  and  structure 
of  the  recirculation  zone. 


(2)  Characterize  the  recirculation  zone  In  terms  of  turbulence  Intensity 


levels  and  velocity  gradients. 

(3)  Investigate  the  fuel  jet/reclrculatlon  zone  Interaction. 

Figure  5a  Illustrates  the  Laser  Anemometry  (LA)  results  for  7.5  m/s, 

^  -  0.05,  Isothermal  (COLD)  flow  and  reacting  (HOT)  flow.  The  formation  of  a 
recirculation  zone  downstream  of  the  centerbody  face  Is  clearly  evident.  The 
fuel  jet  Is  contained  within  the  recirculation  zone  In  the  case  of  Isothermal 
flow.  Reaction  changes  the  flow  structure.  For  the  reacting  case,  the  fuel 
jet  penetrates  the  recirculation  zone.  This  Is  attributed  to  the  higher 
temperature  and  thus  lower  density  of  medium  Into  which  the  fuel  jet  Is 
Issuing.  A  similar  result  was  established  for  the  Isothermal  case  by  doubling 
the  jet  velocity. 

Visually  the  flame  Is  comprised  of  two  regions,  an  upstream  region  with 
the  shape  of  a  seml-elllpsold  attached  to  the  centerbody  face  and  a  downstream 
region  with  the  appearance  of  a  diffusion  flame.  To  the  eye,  the  flame 
appears  steady  and  S3rmmetric  with  the  exception  of  a  rugged  trailing  edge  on 
the  diffusion  flame  and  occasional  oscillations  about  centerline.  Both 
regions  are  luminous  yellow  and  substantial  soot  build-up  occurs  on  the 
centerbody  face.  Indicative  of  a  close  proximity  of  the  reacting  gases  to  the 
solid  boundary. 

The  high-speed  photography  revealed  a  visually  undetected  cyclic  change 
In  the  flame  structure  occurring  at  a  frequency  approximately  110  Hz.  The  no- 
swlrl  flame  oscillates  both  axially  and,  at  a  random  Interval,  radially  about 
the  centerline.  During  each  axial  oscillation  of  the  flame,  a  mass  of 
reacting  fluid  separates  from  the  trailing  edge  of  the  flame,  traveling 
Independently  downstream  In  the  form  of  a  flame  turbule  (Figure  5b).  This 
phenomenon  Is  consistent  with  that  observed  In  the  larger  scale  (AFWAL) 
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Figure  S<  Centerbody  Combustor  Evaluation  (7.5  m/s) 
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system,  combustor  (Roquemore,  et  al.,  1982)  and  Is  attributed  to  vortex 
shedding  off  the  centerbody  face. 

These  observations  established  that,  although  the  no-swlrl  CBC  geometry 
has  value  for  diagnostic  tool  development  (e.g.,  as  employed  at  AFWAL),  not 
all  the  criteria  delineated  above  are  satisfied.  In  particular,  (1)  the 
geometry  does  not  have  swirl,  (2)  the  no-swlrl  flames  exhibit  random  asymmetry 
and  flame  turbule  shedding,  and  (3)  the  flame  Is  attached  to  the  centerbody 
face,  producing  carbon  deposits  and  establishing  an  Intimate  association 
between  the  solid  boundary  condition  and  flame  performance,  a  difficult  If  not 
Impractical  condition  to  model.  This  established  a  point  of  reference  to 
Identify  the  effect  of  swirl  which  was  next  added  to  the  geometry. 

4.1.2  Swirl  CBC 

Swirl  was  added  to  the  CBC  geometry  by  placement  of  45“  straight  blade 
swirl  vanes  upstream  of  the  centerbody  face  (Figure  Ic).  Qualitatively, 
swirling  flows  In  the  CBC  configuration  produce  a  much  more  Intense  and,  to 
the  eye,  stable  flame.  Swirl  tends  to  shorten  the  flame  and  produce  a 
smaller,  but  more  Intense  reaction  zone  with  little  flame  oscillation  apparent 
to  the  eye. 

High-speed  photography  of  the  swirl  flame  revealed  that  the  flame 
undergoes  a  lengthening  and  Shortening  cycle  (on  the  order  of  2:1  In 
magnitude)  at  a  frequency  of  about  100  Hz.  Several  hundred  cycles  were 
observed  (on  the  film)  and,  although  the  period  would  sometimes  vary  ±  50%  In 
duration.  It  was  surprisingly  regular.  The  pictures  also  revealed  that 

(1)  the  swirl  flame  Is  axlsymmetrlc  even  though  It  oscillates  In  length,  and 

(2)  reacting  gases  do  not  separate  from  the  trailing  edge  of  the  flame  during 
these  oscillations,  as  was  the  case  for  the  no-swlrl  flame. 

Figure  6a  shows  the  effect  of  adding  swirl  to  a  7.5  m/s  (bulk  velocity) 
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a)  Axial  and  Radial  Mean  Velocity  Profiles  (7.5  m/s,  COLD,  No  Fuel  Jet] 
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cold  flow  with  no  central  fuel  jet.  Both  the  recirculation  zone  length  and 
structure  experience  significant  change.  With  swirl,  the  recirculation  zone 
extends  beyond  the  20  cm  test  section  with  negative  centerline  velocities 
increasing  monotonically  with  axial  distance  from  the  centerbody  face.  The 
extent  of  recirculation  was  first  shown  by  laser  anemometry  and  later  verified 
via  smoke  traces  to  reach  beyond  the  end  of  the  70  cm  long  duct  and  produce  a 
suction  of  gaseous  material  (combustion  products  and  room  air)  from  the  exit 
plane  to  the  centerbody  face.  In  contrast,  the  no-swirl  case  has  a 
recirculation  zone  contained  within  the  test  section  with  high  negative 
velocities  starting  near  the  centerbody  and  reaching  a  stagnation  point  about 
two  centerbody  diameters  downstream. 

Heat  release  tends  to  suppress  exit  plane  suction  as  illustrated  by 
Figure  6b  which  shows  the  same  swirl  case  with  reaction  (7.3  m/s,  • 

0.05).  However,  an  increase  in  reference  velocity  increases  swirl  strength 
and  can  again  Induce  suction  for  the  reacting,  as  well  as  the  non-reacting 
flows. 

Previous  experimental  analyses  of  complex  flow  combustors  in  which  all 
the  inlet  air  is  swirled  reveal  the  presence  of  indefinitely  long 
recirculation  zones  as  evidenced  by  suction  at  the  exit  plane  of  the  combustor 
(e.g.,  Khalil,  et  al.,  1976,  —  low-intensity  swirl  tube  combustor).  Other 
examples  of  100%  inlet  swirl  configurations  susceptible  to  suction,  include 
hl-lntensity  swirl  tube  combustors  (e.g.,  Appleton  and  Heywood,  1973),  swirl 
dump  combustors  (e.g.,  Buckley,  et  al.,  1980)  and  swirl  expansion  combustors 
(e.g.,  Lilley  and  Rhode,  1982).  Published  modeling  results  do  not,  as  a 
general  rule,  indicate  suction  at  the  exit  plane.  However,  this  may  be  an 
artifact  of  the  modeling  process.  Exit  plane  suction  can  be,  and  is, 
predicted  by  elliptic  flow  codes,  but  the  possibility  is  generally  ruled 
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unrealistic  by  the  modeler  and  excluded  by  setting  the  exit  plane  boundary 
condition  to  preclude  a  reverse  flow. 

A  similar  approach  might  be  employed  experimentally,  that  Is  by 
controlling  the  exit  plane  boundary  conditions  to  suppress  suction.  Flow 
stralghteners ,  for  example,  could  be  positioned  In  the  combustor  to  retain  the 
zone  of  recirculation  upstream  of  the  exit  plane.  In  practical  gas  turbine 
combustors,  however,  this  Is  not  realistic  and  aerodynamic  methods  (l.e., 
radial  jet  Injection  and  dilution  air)  are  employed  to  confine  the  zone  of 
recirculation  to  the  primary  zone  and  thereby  prevent  reverse  flow  from 
extending  to  the  exit  plane. 

Several  development  tests  were  conducted  In  an  attempt  to  eliminate  the 
swirl  Induced  suction  while  maintaining  the  CBC  geometry  as  close  as  possible 
to  the  original  one  described  In  Figure  Ic.  Downstream  swirl  suppression 
devices  such  as  orifice  plates  with  varying  size  holes  were  Installed.  They 
were  found  to  suppress  suction  but  In  doing  so  modify  the  structure  of  the 
recirculation  zone  and  thereby  compromise  flame  stabilization.  Flames  were 
stable  only  at  very  low  equivalence  ratios  and  were  invariably  attached  to  the 
centerbody  face.  Honeycomb  and  other  antl-swlrl  devices  Installed  downstream 
did  not  prevent  suction  but  only  moved  the  point  of  flow  reversal  upstream  to 
the  location  selected  for  the  device. 

4.1.3  Dilute  Swirl  Combustor  (DSC) 

A  swlrl-stablllzed  geometry  was  developed  to  address  the  deficiencies 
observed  with  the  swirl  CBC  geometry  and  satisfy  the  Identified  criteria.  The 
geometries  of  practical  (gas  turbine)  combustors  were  used  as  a  guide  In  these 
development  tests.  The  strong  swirl  required  to  produce  a  stable  flame  must 
be  dissipated  downstream  prior  to  reaching  the  exit  plane  by  the  addition  of 
non-swlrllng  dilution  air.  Adding  this  air  downstream  of  the  recirculation 
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zone  through  radial  wall  jets  would  properly  represent  gas  turbine  combustors 
but  compromise  the  criteria  to  retain  (1)  clean  and  well  prescribed  boundary 
conditions,  and  (2)  clear  optical  access.  Therefore,  an  alternative  approach, 
amenable  to  subsequent  transformation  to  radial  jet  addition,  was  adopted  and 
evaluated  (Figure  Id). 

Swirl  Strength.  The  DSC  was  initially  evaluated  to  determine  the 
ratio  of  swirl  to  dilution  air  (S/D)  required  to  suppress  exit  plane 
suction.  It  was  found  that  if  the  swirl  flow  was  less  than  seventy  (70) 
percent  of  the  total  flow,  no  suction  occurred.  Whereas  the  upper  limit  of 
percent  inlet  swirl  air  is  dictated  by  the  onset  of  suction,  the  lower  limit 
is  bounded  by  the  minimum  necessary  to  produce  swirl-induced  centerline 
recirculation.  A  typical  scenario  is  depicted  in  Figure  7.  At  zero  (0) 
percent  inlet  swirl  air  (Figure  7a,  S/D  =  0),  the  flame  is  similar  to  that  of 
a  no-swirl  CBC,  (i.e.,  the  flame  is  long,  luminous  yellow,  and  attached  to  the 
fuel  tube  face).  At  the  onset  of  swirl-induced  recirculation,  the  flame 
shortens,  becomes  more  intense,  separates  from  the  fuel  tube  face,  and  the 
regions  of  luminous  yellow  soot  are  suppressed  (Figure  7b,  S/D  -  1/4).  At 
fifty  (50)  percent  swirl  air  (Figure  7c,  S/D  -  1/1),  the  flame  is  lOOZ  blue, 
intense  and  compact.  A  50%  inlet  swirl  (i.e.,  a  swlrl/dllution  air  ratio 
(S/D)  of  unity)  was  selected  as  a  baseline  operating  p>olnt.  For  a  swirl-to- 
dllutlon  ratio  of  unity,  the  value  used  in  the  present  study,  the  swirl  number 
(SN)  obtained  by  integrating  across  the  swirl  vanes  is  0.8;  that  obtained  by 
integrating  the  total  inlet  mass  flux  is  0.3. 

Performance.  Radial  and  centerline  profiles  of  axial  mean  and  rms 
velocities  are  shown  in  Figure  8a  for  the  DSC  operating  under  reacting  (HOT) 
flow  at  a  swlrl/dllution  ratio  (S/D)  of  1/1.  Gaseous  propane  is  introduced 
through  a  cone  annular  nozzle  with  an  overall  equivalence  ratio  of  0.1.  A 
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strong,  compact  recirculation  zone  is  produced  and  positive  centerline 
velocities  are  established  well  before  the  exit  plane. 

High  speed  cinematography  of  the  DSC  reveals  a  steady  and  symmetric 
flame.  No  regular  or  periodic  change  in  either  flame  shape  or  luminosity  is 
observed.  This  behavior  Is  attractive  especially  with  respect  to  the 
application  of  time-averaging  modeling. 

Fuel  Injection.  Both  combustor  aerodynamics  and  fuel  injection  play 
an  Important  role  in  dictating  the  flame  intensity,  blow-off  limits,  and  final 
products  of  combustion.  Mellor  (1980),  for  example,  describes  the  use  of 
differential  nozzle  pressure  to  control  the  location  into  which  the  bulk  of 
the  fuel  is  injected  in  the  wake  of  a  stabilizer  disk.  At  low  injection 
pressures,  the  fuel  is  contained  within  the  recirculation  zone  with  a 
concomitant  promotion  of  soot  formation.  At  high  injection  pressures,  the 
fuel  penetrates  the  bulk  recirculation  zone  and  enters  the  shear  layer 
resulting  in  lower  soot  production  and  Increased  flame  stability. 

Tests  of  various  fuel  injection  nozzles  (Figure  8b)  were  conducted 
in  the  present  study  to  identify  the  Influence  of  spray  pattern  in  the  DSC 
flowfleld.  In  general,  it  was  found  that  the  fuel  Injection  pattern  was  a 
principal  determinate  in  the  control  of  blow-off  limits,  flame  shape,  and 
emission  of  soot  particulate. 

Round  fuel  jet  injection  was  first  assessed  for  both  gaseous 
(propane)  and  liquid  (isooctane)  fuels.  For  propane,  jet  diameters  of  1.32  mm 
(0.052-inch)  and  2.03  mm  (0.080-lnch)  ID  were  tested  and,  for  the  liquid 
isooctane,  a  straight  jet  Delavan  air  blast  atomizer  (Model  30609-4)  was 
evaluated.  Results  were  similar  for  both  fuels.  At  lower  equivalence  ratios 
(eg.,  -  0.05),  the  fuels  were  contained  within  the  recirculation  zone 

resulting  in  luminous  soot  production.  As  equivalence  ratio  was  Increased, 
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the  jet  momentum  Increased  and  the  fuel  penetrated  the  recirculation  zone 
resulting  In,  first,  a  diffusion  type  flame  downstream  of  the  recirculation 
zone  and  finally,  blow-off  at  relatively  low  values  of  overall  equivalence 
ratio  ((^  ■0.2). 

Solid  cone  liquid  pressure  atomization  nozzles  (Delavan  WDB  1.5  and 
1.0,  45**)  were  tested  using  Isooctane  at  a  differential  pressure  of  125  psl. 
Injection  of  a  portion  of  the  fuel  Into  the  center  of  the  recirculation  zone 
once  again  resulted  In  luminous  soot  formation  although  less  than  was  observed 
for  straight  jet  Injection. 

Hollow  cone  nozzles  for  both  liquid  (Delavan  WDA  1.5,  80”)  and 
gaseous  fuels  (cone/annular,  and  Delavan  WDA  16,  80”)  were  tested.  For 
propane  and  Isooctane,  these  nozzles  resulted  In  virtually  no  visual  soot 
production  and  extended  the  blow-off  limits  to  equivalence  ratios  greater  than 
1.0.  Laser  anemometry  data  (Figure  8c)  Illustrates  the  zone  of  Injection. 

Cold  flow  Is  selected  for  this  example  to  demonstrate  the  effect  of  the  jet. 
The  effect  of  fuel  Injection  Is  evident  by  comparing  the  radial  profiles  at 
the  1.0  cm  station  for  ■  0.0  (no  jet)  and  ^  >0.1  for  the  cone/annular  gas 
Injection  nozzle.  The  fuel  Is  Injected  Into  a  region  of  high  velocity 
gradients  which  results  In  high  turbulence  levels  and  enhanced  mixing  of  the 
fuel  with  thermally  and  chemically  excited  Intermediates  from  the 
recirculating  core  flow  and  oxygen  from  the  outer  flow. 

4.1.4  Summary 

Criteria  were  established  for  a  common,  laboratory  scale,  complex  flow 
combustor,  and  an  evaluation  was  conducted  of  candidate  geometries.  A  no- 
swlrl  centerbody  combustor  (CBC)  was  first  evaluated.  The  no-swlrl  flames 
exhibit  random  asymmetry  and  vortex  shedding  off  the  trailing  edge  of  the 
flame.  In  addition,  for  many  of  the  no-swlrl  cases,  the  flame  Is  attached  to 
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the  cenCerbody  face,  a  difficult  If  not  Impractical  condition  to  model* 

All  of  the  Inlet  air  to  the  CBC  was  swirled  and  It  was  found  that  this 
often  resulted  In  exit  plane  suction,  a  condition  not  representative  of 
practical  devices,  and  for  which  downstream  boundary  condition  specification 
Is  difficult.  A  development  study  was  conducted  to  establish  a  swirl 
stabilized  geometry  that  controls  both  the  strength  and  extent  of  swirl. 

The  configuration  developed  and  evaluated,  the  dilute  swirl  combustor 
(DSC),  Is  comprised  of  a  combination  swirl,  no-swlrl  Inlet  and  exhibits  clean 
and  well-defined  boundary  conditions  as  well  as  a  detached  symmetric 
recirculation  zone.  Laser  anemometry  measurements  were  made  to  Identify  the 
strength  and  extent  of  the  recirculation  zone  and  baseline  methods  of  fuel 
Injection,  equivalence  ratios  and  swirl  strength  were  Identified. 

Papers  presented  on  this  work  Include  Brum,  Ikloka,  and  Samuelsen  (1981), 
Brum,  Ikloka,  and  Samuelsen  (1982),  and  Brum  and  Samuelsen  (1982a).  A 
manuscript  has  been  submitted  for  publication  (Brum  and  Samuelsen,  1984a). 

More  detailed  Information  Is  provided  In  a  M.S.  Thesis  (Ikloka,  1982)  and  a 
Ph.D.  Dissertation  (Brum,  1983). 

4.2  Element  B;  Physical  and  Chemical  Processes  In  Complex  Flows 

Attention  under  this  element  was  directed  In  the  present  program  to 
(1)  the  acquisition  and  Interpretation  of  detailed  flowfleld  data  In  the 
premixed  Opposed  Jet  Combustor  (OJC)  and  the  non-premlxed  Dilute  Swirl 
Combustor  (DSC),  and  (2)  the  performance  of  elliptic  flow  modeling  In 
predicting  the  acquired  experimental  data.  Experiments  were  conducted  for  the 
following  conditions:  COLD  (non-reacting),  HEATED  (non-reacting),  and  HOT 
(reacting).  The  justification  of  conducting  non-reacting  flow  experiments 
(COLD  and  HEATED)  was  to  develop  an  experimental  data  base  for  testing 
elliptic  flow  models  In  the  absence  of  reaction.  Such  tests  are  both  useful 
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and  necessary  to  detect  shortcomings  of  the  models  prior  to  tests  for  the 
cases  In  which  the  complication  of  reaction  Is  added. 

4.2.1  Opposed  Jet  Combustor  (COLD  Experiments) 

Two  COLO)  non-reacting  opposed  Jet  combustor  experiments  were  conducted 
under  the  present  program.  One  was  conducted  In  the  Laboratory  Combustor 
Facility  (LCF)  and  the  second  In  the  Hydrodynamic  Facility  (HOF). 

LCF  Tests 

The  opposed  Jet  module,  depicted  In  Figure  la,  was  operated  In  a 
nonreacting  (Isothermal)  mode  wherein  the  main  stream  flow  was  air  and  the  Jet 
stream  flow  was  carbon  monoxide.  The  baseline  main  stream  and  Jet  stream 
velocities  were  15.24  m/sec  (50  ft/sec)  and  140  m/sec  (460  ft/sec) 
respectively.  The  Reynolds  and  Mach  numbers  of  the  main  flow  were  50,000  and 
0.05  respectively,  and  11,000  and  0.4  respectively  for  the  Jet  flow.  Reynolds 
numbers  were  sufficiently  high  to  Insure  turbulent  flow  In  both  the  pipe  and 
Jet,  and  both  flow  sources  had  entrance  lengths  sufficient  to  provide  fully 
developed  pipe  flow  conditions.  Mach  numbers  were  sufficiently  low  to  avoid 
compressibility  effects.  This  operating  condition  produced  a  recirculation 
zone  which  extended  approximately  4  cm  (0.78  combustor  diameters)  upstream  of 
the  Jet  exit  plane.  A  second  operating  condition  was  established  In  which  the 
velocity  of  the  main  stream  was  reduced  by  50  percent  to  7.62  m/sec  (25 
ft/sec).  This  operating  condition  resulted  In  a  substantially  larger  upstream 
extension  of  the  recirculation  zone  (6.3  cm,  or  1.2  combustion  diameters). 

The  objective  of  the  experimental  measurements  was  to  provide  a  basis  for 
Judgment  of  modeling  results  as  well  as  to  assist  In  the  specification  of  the 
flow  Inlet  boundary  conditions.  The  key  diagnostics  Included: 
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(1)  a  laser  anemometer  for  the  measurement  of  Instantaneous  velocity  and 
the  subsequent  determination  of  time-averaged  velocity  and 
turbulence  intensity, 

(2)  a  pitot  probe  for  the  measurement  of  mean  velocity  to  provide  a 
check  for  and  extend  the  laser  anemometer  measurements,  and 

(3)  a  nondlsperslve  infrared  analyzer  and  extractive  probe  sampling 
system  for  the  measurement  of  the  diffusion  of  the  tracer  (carbon 
monoxide)  throughout  the  flowfleld. 

The  specific  data  developed  were  axial  velocity  profiles  at  ten  stations 
within  the  combustor,  contours  of  constant  turbulence  Intensity,  and  contours 
of  constant  mass  concentration  of  the  carbon  monoxide  tracer  species. 

Two  sets  of  experimental  data  were  obtained:  (1)  momentum  transport 
measurements  using  laser  anemometry  and  pitot  probe  techniques,  and  (2)  mass 
transport  measurements  defining  the  diffusion  of  the  CO  tracer  species  using 
probe  sampling  and  nondlsperslve  infrared  detector  measurements. 

Laser  anemometer  measurements  were  restricted  to  that  region  forward  of 
the  jet  exit  plane  because  of  the  Interference  of  the  jet  body  with  the  laser 
beams.  The  key  results  of  the  laser  anemometer  measurements  obtained  were: 

(1)  definition  of  the  combustor  Inlet  velocity  and  turbulence  Intensity 
profiles , 

(2)  definition  of  the  flow  development  In  the  main  combustor  tube,  and 

(3)  definition  of  the  velocity  and  turbulence  Intensity  In  the  region 
surrounding  and  In  the  outer  portion  of  the  recirculation  zone. 

Figure  9a  shows  the  Inlet  velocity  profile  at  the  combustor  Inlet  for  the 
baseline  15.24  m/sec  main  flow  test  condition.  Shown  are  two  full  diameter 
profiles  folded  about  the  axis  of  symmetry  to  display  both  profile  symmetry 
and  repeatability. 
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RADIUS 


a)  Inlet  Velocity  Profiles 


Mass  diffusion  was  measured  within  the  recirculation  zone  and  the  wake. 

Proper  design  of  sample  probe  configurations  allowed  mapping  those  regions  of 

the  flowfleld  where  the  concentration  ranged  from  0.5  to  8  percent  (Figure 

9b).  Numerical  results  were  obtained  for  the  (1)  velocity  components  (u  and 

v),  (2)  turbulence  kinetic. energy  (k),  (3)  turbulence  kinetic  energy 

dissipation  rate  (e)  ,  (4)  effective  viscosity  (y  and  (5)  tracer  species 

et  t 

(carbon  monoxide)  concentration  (CO). 

The  sensitivity  of  the  numerical  solution  to  various  features  of  the 
analytical  models  (e.g.  transport  models,  model  parameters,  and  boundary 
conditions)  were  Investigated.  Specifically,  the  following  model  requirements 
were  addressed: 


Case 

Subject 

(1) 

turbulence  model  specification. 

(2) 

density  of  the  numerical  grid 

(3) 

turbulence  kinetic  energy  dissipation  rate  on  jet  body  face, 

(4) 

Inlet  turbulence  kinetic  energy  for  the  main 

and  jet  flows. 

(5) 

Inlet  turbulence  energy  dissipation  rate  for 

the  main  and 

jet  flows. 

(6) 

Inlet  velocity  profiles  for  the  main  and  jet 

flows, 

(7) 

variations  In  main  and  jet  flow  rates,  and 

(8) 

turbulence  model  parameters  (Cp  C2,  C  ,  C^). 

Velocity  and  species  concentration  predictions  could  be  directly  compared 
to  the  experimental  measurements.  As  a  result,  these  comparisons  were  used  as 
the  primary  means  of  evaluation  of  the  analytical  results.  The  findings  are 
summarized  In  Table  I  and  discussed  below  by  case. 


39 


Case  1 .  The  algebraic  viscosity  model  did  not  exhibit  the 
capability  to  adequately  predict  the  complex  flowfleld  of  the  combustor 
(Figure  10).  The  resultant  constant  turbulent  viscosity  could  not 
simultaneously  produce  the  trends  observed  experimentally  of  jet  penetration 
and  mass  diffusion  in  the  wake.  Solutions  using  this  model,  however,  were 
valuable  in  pointing  out  the  magnitude  of  the  turbulent  viscosity  needed  to 
produce  the  diffusion  of  momentum  and  mass  observed  experimentally  and 
facilitating  a  first  solution  with  a  minimum  of  independent  variables. 

The  two-equation  k-e  model  did  produce  a  distribution  of  turbulence 
within  the  flowfleld,  and  displayed  the  potential  to  more  satisfactorily 
predict  the  combustor  flowfleld.  However,  certain  deficiencies  were  apparent 
in  the  ability  of  the  model  to  predict  experimental  trends.  For  example: 

(1)  The  velocity  profiles  (Figure  10a)  show  that  the  predicted  momentum 
exchange  between  the  main  stream  and  the  recirculation  zone  is 
insufficient.  (This  result  is  evident  from  inspection  of  the 
velocity  profile  at  the  station  X/Rq  ■  1.0). 

(2)  The  predicted  velocity  profiles  adjacent  to  the  outer  wall  of  the 
jet  display  a  momentum  deficit. 

(3)  The  predicted  concentration  gradients  (Figure  10b)  in  the  vicinity 
of  the  stagnation  region  are  significantly  steeper  than  those 
measured. 

(4)  The  predicted  radial  diffusion  of  mass  in  both  the  recirculation 
zone  and  the  wake  was  limited. 


TABLE  I 


-  Summary  of  Opposed  Jet  COLD  Modeling  Study  - 


Caae 

Sensitivity  Investigated 

Results 

1. 

Turbulence  model 

algebraic  model 

U  •  constant,  could  not  simultaneously  match  penetra- 
and  radial  diffusion  of  Jet 

k-e  model 

limited  radial  diffusion  of  momentum  and  mass,  high 
concentration  gradients  predicted 

2. 

Numerical  grid 

18x31  grid 

4  horizontal  gridlines  in  Jet,  9  vertical  gridlines 

In  recirculation  zone,  provided  a  realistic  solution 

25x49  grid 

14  horizontal  gridlines  In  Jet,  22  vertical  gridlines 
In  recirculation  zone,  provided  Improved  solution 

3. 

Recirculation  zone 
wall  boundary 
conditions 

0.1  £  (nom) — 
w 

1.0  e  (nora) 
w 

.  high  sensitivity,  solution  develops  unrealistic 
character  for  specification  of  off  nominal  condition 

4, 

Inlet  turbulence 
Iclnetlc  energy 

k(nom)  ”  .0054  U* 
l,0k(nom)~6.0k(noiJi) 

no  significant  sensitivity  noted 

5. 

Inlet  turbulence 
kinetic  energy 
dissipation  rate 

e  »  k^/^/1 

2 (nom)  •  O.IOR 
.051 (nom) --1.01 (nom) 

no  significant  sensitivity  noted 

6. 

Inlet  velocity 
profiles 

square 

realistic  solution 

turbulent  pipe 
profiles 

no  significant  change 

7. 

Main  and  Jet  flow 
rate  variations 

0.9U^(nom)  and 

no  significant  change 

I.IU  (nom) 
m 

0.5^  (nora)  — 
l.OiJn  (nom) 

significant  increase  in  extent  of  recirculation  zone 
as  main  stream  velocity  decreased 

l.OC  (nom) — 

3. 

Turbulence  model 
parameters 

a.  0.5  Cj^(nom) 

b.  1.5  C2(nora) 

c.  0.5  Cjj(nom 

d.  1.5  C  (nora) 

u 

Jet  penetration  decreased,  radial  diffusion  increas¬ 
ed,  and  stagnation  zone  gradients  decrease  for  all 

cases 

f.  0.875  C.(nom) 

most  improved  solution 

a)  Momentum  Transport 


U/0,  •  MEASURED  VALUE  ^  PREDICTED  VALUE 


18  X  31  GRID,  ALGEBRAIC  TURBULENCE  MODEL,  SQUARE  INLET  VELOCITY  PROFILES 


AXIAL  STATION  X/R„ 


b)  Mass  Transport 


CO  TRACER  CONCENTRATION  PROFILES 


Figure  10.  Opposed  Jet  COLD  Combustor  Flowfield  Data: 
Sensitivity  to  Turbulence  Model 
(15.24  m/s) 
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These  results  indicate  that  the  magnitude  of  the  predicted  turbulent 
viscosities  is  underestricted  and  that  the  spatial  distribution  of  viscosity 
is  not  adequately  predicted. 

Case  2.  Refinement  of  the  computational  grid  to  provide  greater 
detail  within  the  jet,  near  the  Jet  exit,  and  within  the  recirculation  zone 
resulted  in  a  significant  change  in  the  character  of  the  resultant  solution. 
The  changes  reflect  an  Improved  match  of  the  predicted  and  experimental 
results  (Figure  11),  and  result  from  the  reduction  of  numerical  error  in  the 
application  of  the  finite  difference  procedure. 

Case  3.  Previous  studies  have  noted  that  in  confined  flows,  where 
the  recirculation  zone  may  contact  solid  boundaries,  the  specification  of 
near-wall  turbulence  energy  dissipation  rate  may  significantly  influence  the 
mixing  characteristics  (Peck  and  Samuelsen,  1977).  This  proposal  was  tested 
in  the  present  investigation  both  experimentally  and  numerically  by  reducing 
the  outer  diameter  of  the  jet  tube  from  6.4  to  1.6  mm. 

Experimentally  measured  profiles  of  turbulence  Intensity  mass  show 
less  turbulence  and  less  mass  mixing  for  the  case  where  the  small  outside 
diameter  jet  was  used.  The  differences,  however,  are  quite  small  and  do  not 
support  the  proposal  that  the  presence  of  the  jet  body  face  in  the  vicinity  of 
the  recirculation  zone  has  a  large  Influence  on  the  structure  of  the  resultant 
flowfleld. 

In  contrast,  the  predictions  were  found  to  be  quite  sensitive  to  the 
specification  of  the  energy  dissipation  rate  on  the  upstream  facing  step  of 
the  jet  tube.  The  use  of  other  than  the  nominal  value  of  the  energy 
dissipation  rate  on  solid  boundaries  does  not  appear  justified  and  may  in  fact 
mask  fundamental  shortcomings  of  the  model. 

Case  4  and  5.  The  predictions  were  found  to  be  not  sensitive  to  the 
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a)  Momentum  Transport 


U/U, 
0  1.0 


•MEASURED  VALUE 

0  1.0  n  I 


— -  PREDICTED  VALUE 


18  X  31  GRID,  2-EQUATrOM  TURBtn.ENCE  MODEL.  SQUARE  INLET  VELOCITY  PROFILES 


25  X  49  GRID,  2-EqUATION  TURBULENCE  MODEL,  SQUARE  INLET  VELOCITY  PROFILES 

AXIAL  STATION  X/R„ 


b)  Mass  Transport 


CO  TRACER  CONCENTRATKM  PRORLES 
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Figure  11.  Opposed  Jet  Combustor  COLD  Flowfield  Data: 
Sensitivity  to  Numerical  Grid 
(15.24  m/s) 
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specification  of  inlet  turbulence  kinetic  energy  and  inlet  turbulence  energy 
dissipation  rate.  Comparable  magnitudes  of  turbulence  Intensity  were 
developed  in  the  zone  of  primary  activity  (i.e.  recirculation  zone)  for  a 
sizeable  range  of  values  for  both  parameters.  The  reason  for  the 
insensitivity  observed  is  attributed  to  the  separation  of  the  recirculation 
zone  from  the  flow  entrance  planes  in  terms  of  pipe  diameters  and/or  jet 
diameters.  As  a  result,  the  viscosity  model  dominates  the  flowfleld 
development  rather  than  the  flow  inlet  boundary  conditions.  Were  the 
recirculation  zone  located  in  close  proximity  to  the  flow  Inlets ,  the  solution 
would  be  expected  to  be  sensitive  to  the  inlet  flow  boundary  conditions. 

Case  6.  To  add  to  the  realism  of  the  boundary  condition 
specification,  velocity  profiles  were  input  for  both  the  main  and  jet  flow 
inlets.  The  measured  velocity  profile  (Figure  9a)  was  input  for  the  main  flow 
and  a  one-seventh  power  law  turbulent  pipe  flow  was  input  for  the  jet  flow. 

No  significant  change  in  the  resultant  solution  occurred  as  a  result  of  this 
refinement.  The  reason  for  the  insensitivity  to  this  refinement  is  identical 
to  that  given  for  the  insensitivity  to  inlet  k  and  e  specification.  That  is, 
the  primary  zone  of  activity  (i.e.,  the  recirculation  zone)  is  sufficiently 
removed  from  the  flow  entrance  prolanes  in  terms  of  pipe  diameters  and/or  jet 
diameters.  Were  this  not  the  case,  the  solution  would  be  expected  to  be 
sensitive  to  the  specification  of  the  inlet  momentum  distribution. 

Case  7.  The  sensitivity  of  the  flowfleld  to  both  small  and  large 
changes  in  input  flow  momentum  was  investigated.  Small  changes  in  input  flow 
momentum  were  first  investigated  to  examine  the  potential  Influence  of 
experimental  flow  rate  errors  on  the  resultant  solution.  Insignificant  change 
in  the  solution  was  observed  for  a  simultaneous  +10  percent  change  in. the  main 
flow  rate  and  -10  percent  change  in  the  jet  flow  rate.  These  perturbations 
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significantly  exceed  any  errors  present  in  the  experimental  procedure  and  show 
an  insensitivity  of  the  experimental  results  to  potential  flow  rate  errors. 

A  significant  change  in  the  flowfleld  character  resulted  when  the 
main  flow  velocity  was  reduced  to  50  percent  from  15.24  m/sec  to  7.62  m/sec. 
The  axial  extent  of  the  recirculation  zone  approximately  doubled  as  a  result 
of  the  perturbation.  The  same  general  disagreement  between  prediction  and 
experiment  was  noted  for  the  7.62  m/sec  flow  as  for  the  15.24  m/sec  flow. 
Specifically,  the  extent  of  the  radial  diffusion  of  momentum  and  mass  was  less 
for  the  prediction  and  concentration  gradients  in  the  vicinity  of  the 
stagnation  point  were  steeper.  These  results  demonstrate  that  the  disparities 
between  experiment  and  prediction  were  consistent  for  a  wide  range  of  flow 
conditions. 

Case  8.  The  predictions  developed  thus  far  all  indicate  common 
points  of  disagreement  with  the  experimental  data.  The  predictions  for  both 
the  diffusion  of  momentum  and  the  diffusion  of  mass  point  to  an 
underestimation  of  the  turbulent  viscosity  within  the  recirculation  zone. 

This  is  indicated  by  the  greater  predicted  axial  momentum  of  the  flow  in  that 
portion  of  the  recirculation  zone  where  velocities  are  small  and  by  the 
limiting  of  the  radial  diffusion  of  mass  within  this  same  region. 

The  key  to  Improved  solutions  appears  to  lie  in  an  Increase  in 
turbulence  within  the  recirculation  zone.  It  has  been  pointed  out  that  the 
turbulence  model  itself  can  be  tuned  to  accomplish  such  effects  (Launder  and 
Morse,  1977).  In  the  present  study,  the  sensitivity  of  the  resultant  solution 
to  pertubations  of  various  coefficients  of  the  viscosity  model  was  examined 
(Table  I).  It  was  found  that  significant  improvements  in  the  solution 
resulted  from  decreasing  the  coefficient  of  the  source  term  in  the  energy 
dissipation  equation  (Figure  12).  The  result,  less  turbulence  energy 
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dissipation,  produced  greater  turbulence  levels  and  the  enhanced  mixing  of 
mass  and  momentum  throughout  the  entire  flowfleld. 

Summary .  A  turbulence  model  which  accounts  for  the  spatial 
generation  of  and  dissipation  of  turbulence  Is  essential  to  obtaining  a 
realistic  solution.  The  two-equation  k-e  model  shows  limited  promise  In 
describing  the  complex  flow  adopted  for  the  present  Investigation.  The  major 
shortcoming  Identified  for  the  case  of  nonreacting  (COLD)  flow  Is  In  the 
underprediction  of  turbulent  viscosity  within  the  recirculation  zone  and  the 
wake.  Detailed  Information  Is  available  In  a  Ph.D.  Dissertation  (Wuerer, 
1978).  One  paper  was  presented  on  this  work  (Wuerer  and  Samuelsen,  1979),  and 
the  Information  Is  being  Incorporated  Into  a  manuscript  planned  for  submission 
to  Experiments  In  Fluids  (Samuelsen,  et  al,  1984). 

HDF  Tests 

The  results  of  the  OJC  tests  In  the  LCF  described  above  and  those 
described  below  In  Section  4.2.2  made  It  clear  that  (1)  a  substantial  advance 
Is  necessary  In  turbulence  modeling  for  accurate  predictions  of  complex  flows, 
and  (2)  a  need  exists  for  spatially-resolved  measurements  In  complex  flows  for 
conditions  ranging  In  complexity  from  non-reacting  (COLD  and  HEATED)  to 
reacting  (HOT).  Such  data  were  required  to  provide  the  foundation  upon  which 
numerical  modeling  can  be  developed  and  tested,  and  to  provide  the  Insight 
Into  the  physical  and  chemical  processes  of  transport  and  reaction  as  a  guide 
to  both  the  evolution  of  modeling  and  the  understanding  of  these  type  of 
flows. 

A  second  vehicle  selected.  In  addition  to  the  LCF,  for  acquisition  of 
these  data  was  the  Hydrodynamic  Facility  (HDF).  Data  were  obatlned  with  the 
opposed  jet  configuration  In  the  Hydrodynamic  Facility  (HDF)  for  two 
reasons.  First,  such  data  provide  a  data  base  amenable  to  Isothermal 
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modeling.  Second,  water  flows  are  more  amenable  than  gaseous  flows  to  flow 
visualization,  the  latter  of  which  can  be  used  as  a  check  of  the  experimental 
data  as  well  as  to  provide  Insight  Into  both  the  mixing  and  dynamics  of  the 
flowfleld. 

The  opposed  jet  flowfleld  was  Investigated  In  the  HDF  for  nominal 
approach  velocities  of  0.32  m/s  and  0.64  m/s  at  a  constant  jet  velocity  of  5.8 
m/s.  As  shown  In  Table  II,  these  velocities  were  selected  to  yield  the 
Identical  Reynolds  numbers  of  both  the  mainstream  and  jet  examined  In  the 
LCF.  The  variation  In  mainstream  approach  vellclty  provides  a  means  of 
examining  recirculation  zone  strength  and  also  demonstrates  the  Influence  on 
experimental  and  analytical  trends  of  large  changes  In  differential  flowrate 
between  the  main  and  jet  streams. 

Representative  experimental  and  predicted  profiles  of  the  axial  component 
of  mean  velocity  are  presented  In  Figure  13  for  both  mainstream  flowrates. 
Figure  13a  presents  centerline  velocity  profiles,  while  Figure  13b  provides  a 
radial  profile  at  an  axial  location  In  the  recirculation  zone  near  the 
stagnation  point. 

The  trends  exhibited  In  the  analytically  determined  mean  velocity 
profiles  using  a  18  x  31  grid  and  plug  velocity  profile  In  the  jet  agree 
qualitatively  with  those  established  experimentally.  The  axial  penetration  of 
the  opposing  jet  as  determined  numerically,  however,  extends  farther  upstream 
than  Indicated  experimentally;  this  effect  Is  observed  for  both  the  low 
Reynolds  number  (0.32  m/s  average  velocity)  and  the  high  Reynolds  number  (0.64 
m/s  average  velocity)  malnflow  cases.  The  jet  expansion  also  Is  overestimated 
by  the  analytical  model  In  both  cases. 


TABLE  II 


-  HDF  Test  CcndlClons  - 


MAINSTREAM 

JETSTREAM 

HDF 

HDF 

OJC 

HDF 

HDF 

OJC 

Velocity 

Flowrate 

Velocity 

Re 

Velocity 

Flowrate 

Velocity 

Re 

(m/s) 

(  i/m) 

(m/s) 

(m/s) 

(gpm) 

(m/s) 

0.32 

87 

7.5 

25,000 

5.8 

0.28 

135 

11,500 

0.64 

174 

15 

50,000 

5.8 

0.28 

136 

11,500 

50 


b)  Radial  Profiles 
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The  experimentally  determined  jet  penetration  distance  Is  8.4  cm  (1.1 
pipe  diameters)  for  the  0.32  m/s  mainstream  approach  velocity;  In  contrast, 
the  numerical  code  predicts  a  stagnation  point  at  9.2  cm  (1.2  pipe 
diameters).  LA  measurements  for  the  high  Reynolds  number  case  Indicate  a  jet 
penetration  of  4.2  cm  (0.55  pipe  diameters)  In  the  axial  direction;  this 
result  compares  with  a  numerically  predicted  value  of  5.4  cm  (0.71  pipe 
diameters).  In  Table  III,  these  data  are  compared  to  those  obtained  In  the 
LCF.  The  overestimate  of  jet  effects  In  the  HDF  Is  consistent  with  those 
observed  In  the  LCF  for  a  course  numerical  grid  and  plug  velocity  profile  in 
the  jet. 

Finally,  flow  visualization  (Figure  13c)  using  a  dye  In  the  jet  serves  as 
a  check  of  the  LA  data  and  graphically  Illustrates  the  jet  penetration  and 
size  of  the  recirculation  zone. 

Additional  details  are  provided  In  a  M.S.  Thesis  (Acroumanls,  1980)  and  a 
M.S.  Project  Report  (Grudovlch,  1983). 

4.2.2  Opposed  Jet  Combustor  (HEATED  Experiments) 

In  the  work  described  above,  tests  for  the  description  of  mass  and 
momentum  transport  were  conducted  In  the  absence  of  reaction  and  heat 
release.  Heat  transport  was  assessed  by  heating  the  jet  In  the  LCF  and 
comparing  predicted  values  of  mean  and  rms  temperature  with  experimentally 
measured  values.  The  opposed  jet  was  operated  for  this  study  at  rounded 
values  of  the  mainstream  velocity,  7.50  and  15.0  m/s.  Instead  of  the  7.62  and 
15.24  m/s  reference  velocities  used  above. 

The  temperature  signals  were  obtained  by  means  of  0.125  mm  diameter  glass 
coated  thermistors  and  a  1.25  mid  diameter  resistance  ("cold  wire") 
thermometer.  The  cold  wires  were  platinum  and  0.48  mm  In  length  for  ~ 

7.5  m/s,  and  platlnum-10  percent  rhodium  and  0.66  mm  In  length  for  Uq  ^  ~ 


TABLE  III 


FACILITY 


HDF 


LCF 


-  Location  of  OJC  Stagnation  Point  - 


NUMBER  of  pipe  diameters  from  Jet  exit 


TECHNIQUE 


Re  -  25,000 


Re  -  50,000 


LA 

1.1 

0.55 

Prediction 

1.2 

0.71 

Flow  Visualization 

0.92 

0.66 

LA 

1.2 

0.60 

Prediction 


1.3 


0.70 


15.0  m/s.  The  cold  wires  were  operated  with  a  root  mean  square  current  of  255 
microampers. 

The  radial  profiles  of  0/0  and  T’/0  at  x/d  *  2.56  for  the  7.5  m/s 

max  max 

flow  are  shown  in  Figure  14a  where  0  is  the  time-mean  temperature  difference 

(T  -  Tj)  between  the  local  and  inlet  temperature.  (Because  of  interference 

between  the  probe  and  jet  body,  measured  values  were  not  obtained  at  radial 

locations  of  r/R  less  than  0.15.)  The  predictions  depict  a  peak  at  r/R  of  0.1 

in  the  T'/0  distribution;  this  is  associated  with  the  steep  gradient  in 
max 

0/0  at  that  radial  location.  In  the  outer  region,  the  predictions  are  in 
max 

fair  agreement  with  measurements. 

Figure  14b  shows  the  radial  profiles  of  0/0  and  T'0  at  four  axial 
®  max  max 

locations  for  the  15  m/s  flow.  The  axial  and  radial  thermal  extent  of  the 

heated  jet  is  shown  on  both  figures.  It  is  Interesting  to  note  that,  as 

expected,  the  radial  extent  of  the  temperature  fluctuations  is  larger  than 

that  of  the  mean  temperature  at  all  the  four  stations. 

The  mean  temperature  is  well  predicted  at  stations  farther  from  the  jet 

exit  and  overpredicted  near  it.  This  may  be  attributed  to  overpredicted 

turbulent  diffusion  coefficients  which  could  be  due  to  underpredicted  values 

of  the  rate  of  dissipation  of  kinetic  energy,  e  ,  or  due  to  the  use  of 

constant  o  . 

H 

At  the  station  x/d  >  2.95,  the  predicted  0/ 0  is  in  excellent  agreement 

max 

with  the  experimental  data.  However,  at  the  same  station,  the  temperature 
fluctuations  are  underpredicted. 

It  should  be  mentioned  here  that  the  transport  equation  for  the  mean 
square  fluctuation  of  a  scalar  quantity  (LaRue,  Deaton,  and  Gibson,  1975}  has 
been  validated  for  turbulent  free  jets  (Lockwood  and  Naguib,  1975)  and  for 
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confined  turbulent  recirculating  flows  (Elghobashi,  Pun,  and  Spalding, 

1977).  The  same  equation  (with  both  the  values  of  C^2  2.0  and  1.4)  did  not 

accurately  predict  the  measured  values  of  rms  temperature  fluctuations  under 
the  experimental  conditions  of  this  study.  This  stresses  the  need  to  solve  a 
transport  equation  for  the  dissipation  of  the  temperature  fluctuation.  These 
data  also  reinforce  the  need  for  a  closer  examination  of  the  e  equation  and 
the  assumption  of  the  constant  turbulent  Prandtl  and  Schmidt  numbers.  One 
such  example  is  Hanjallc,  Launder,  and  Schlestel,  1979. 

Detailed  information  of  this  study  is  available  is  publication 
(Elghobashi,  et  al,  1981). 

4.2.3  Opposed  Jet  Combustor  (HOT  Experiment) 

The  objective  of  the  present  study  was  to  provide  a  general 
characterization  of  the  chemical  and  flame  structure  of  the  OJC  (Figure  la) 
for  a  range  of  parametric  variations  of  the  four  primary  controlling 
variables:  Main  and  jet  stream  equivalence  ratios,  main  and  jet  stream 
velocities.  Exit  plane  and  detailed  flowfield  profiles  are  presented  and 
analyzed  for  NO^^,  carbon  monoxide  (CO),  total  hydrocarbons  (HC),  and 
temperature.  The  goals  are  to  provide  (1)  insight  into  the  performance  of  a 
reverse  jet,  aerodynamic  flameholder,  (2)  guidance  for  practical  applications 
of  aerodynamic  flameholding,  and  (3)  a  data  base  for  future  code  testing. 

For  detailed  flowfield  maps,  radial  traverses  were  taken  at  twelve  axial 
locations,  with  the  distance  between  axial  locations  ranging  from  2.54  mm 
(0.1-inch)  in  regions  of  steep  gradients  (e.g.,  the  nose  of  the  recirculation 
zone)  to  up  to  24.13  mm  (0.95-lnch)  in  regions  of  near  zero  gradients  (e.g., 
downstream  of  the  jet  tube  exit).  Radial  locations  were  at  3.05  mm 
(0.12-inch)  increments  between  the  Jet  tube  wall  and  the  chamber  wall. 

The  test  conditions  are  tabulated  in  Table  IV.  The  following  conditions 


58 


e 


TABLE  IV 

-  OJC  HOT  Test  Conditions  - 


Uj(m/s) 

cn 

Parametric  variation 

7.5 

135 

1.0 

1.0 

BASE  CASE 

7.5 

135 

1.2 

0.8, 

1.0, 

1.2 

1.0 

0.8 

0.8, 

0.8, 

1.0, 

1.0, 

1.2 

1.2 

0.6 

0.8, 

1.0, 

1.2, 

1.4,  1.6 

1.0  1.0 

1.0  1.0 


7.5 

15.0 
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135 


represented  the  BASE  CASE: 

•  Main  Stream  Reference  Velocity  U  ■  7.5  m/s 

m 

•  Jet  Stream  Reference  Velocity  *  135  m/s 

•  Main  Stream  Equivalence  Ratio  1.0 

•  Jet  Stream  Equivalence  Ratio  >■  1.0 

Results  are  presented  for  the  BASE  CASE.  Second,  parametric  variations 

are  presented  In  the  order  listed  in  Table  IV.  Each  of  the  four  controlling 

variables  (*  ,  6.,  U , ,  U  )  were  evaluated.  Note  (Table  IV)  that  the 
m  j  j  m 

mainstream  equivalence  ratios  (i^^)  were  biased  to  fuel  lean  mixture  ratios 
because  of  the  Interest  In  lean  mainstream  emission  performance.  The  Jet 
equivalence  ratios  ( ()ij )  were  biased  to  the  fuel  rich  mixture  ratios  to  extend 
the  lean  limit  of  the  main  stream  mixture.  The  range  of  main  stream  (U^)  and 
jet  velocities  (Uj)  allowed  an  examination  of  the  effects  of  recirculation 
zone  size  and  stoichiometry  on  flame  structure  and  pollutant  emission. 

Finally,  emission  Indexes  are  presented  to  summarize  the  emission 
behavior  of  the  combustor  at  the  conditions  considered,  and  to  provide  a 
practical  perspective  to  the  utility  of  opposed  Jet  flameholding. 

Base  Case.  The  detailed  flowfleld  maps  and  exit  plane  profiles  for 
the  base  case  are  presented  In  Figure  15a.  Two  distinct  regions  can  be 
deduced  from  the  results  (Figure  15b).  One  Is  the  recirculation  zone,  which 
Is  a  zone  of  strong  backmlxlng  driven  by  the  Jet  flow.  The  other  Is  a 
radially  propagating  reaction  In  the  wake  of  the  recirculation  zone.  For 
example,  the  oxidation  of  hydrocarbons  and  formation  of  carbon  monoxide  (CO) 
occurs  within  the  recirculation  zone  where  there  Is  Intense  mixing  of 
reactants  with  hot  products,  and  along  the  radially  propagating  wake  reaction 
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b)  Heat  Release  and  Pollutant  Formation  Scenario 


HC  oxidized 
CO  formed 


HC  oxidized 


c)  Parametric  Study:  Effect  of  U  and  U.  (Flame  Shape) 

_ ; _ !s _ j _ 


Base  case 


M| 

1 , - - 

Elevated  mainstream  velocity  {V,») 


i 


Figure  15.  (continued) 


g)  Parametric  Study:  Effect  of  S  and  6.  (N0_) 

a  j  * 


h)  Parametric  Study:  Effect  of  ^  and  ,  (NO/NO  ) 


Figure  IS.  (continued) 
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NO/NOx 


Chamber  Wall 


front.  Within  the  wake,  temperature,  oxygen,  and  residence  time  are 
sufficient  to  ensure  nearly  complete  HC  consumption,  and  to  initiate  the 
oxidation  of  CO  to  carbon  dioxide  (CO2}  as  demonstrated  by  the  decrease  in  CO 
concentration  adjacent  to  the  jet  tube  proceeding  downstream  toward  the  exit 
plane.  Proceeding  toward  the  combustor  wall,  the  concentrations  of  HC  and  O2 
approach  those  of  the  reactants.  As  a  result,  the  source  of  the  hydrocarbons 
emitted  at  the  exit  plane  is  the  area  outside  of  the  wake. 

Oxides  of  nitrogen  (NO^^)  are  formed  thermally  in  both  the 
recirculation  zone  and  wake  as  a  result  of  elevated  temperatures,  sufficient 
residence  time,  and  available  oxygen.  Area-averaged  concentrations  calculated 
at  both  the  “jet  exit  plane”  and  combustor  “exit  plane"  (Figure  la)  Indicate 
that  75  percent  of  the  total  NO^  emitted  is  formed  in  the  recirculation  zone 
for  this  base  condition. 

The  exit  plane  profiles  show  the  general  structure  of  the  wake. 
Within  the  wake  and  proceeding  from  the  jet  tube  to  the  combustor  wall,  HC  and 
oxygen  concentrations  and  temperature  remain  relatively  constant,  while  CO 
concentrations  Increase  slowly  and  concentrations  decrease.  At 
approximately  r/R  «  0.55,  the  concentrations  of  HC,  CO2,  and  O2  change 
sharply.  Oxygen  and  HC  rise,  CO2  drops  and  CO  peaks.  Eventually,  the  HC  and 
O2  rise  to  the  reactant  concentrations. 

Finally,  it  is  noteworthy  that  the  NO/NO^^  ratio  drops  abruptly  at 
the  flame  front.  This  is  attributed  to  the  rapid  mixing  of  hot  products  and 
cold  reactants  at  the  flame  front  which  produce  radical  relaxation  reactions 
and  associated  populations  of  hydroperoxy  radicals  (HO2)  sufficient  to  oxidize 
NO  to  NO2.  Unfortunately,  these  events  can  be  Influenced  by  the  probe,  and 
the  extent  to  which  the  measured  levels  of  NO2  are . real  or  artifacts  of  the 
probe  remains  unanswered.  However,  an  evaluation  (Chen,  et  al.,  1979)  of 
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similar  observations  In  a  premixed  combustor  (Oven,  Gouldln,  and  McLean,  1979) 
concluded  that,  although  measurements  within  high  temperature  reaction  zones 
(e.g.,  within  the  recirculation  zone  and  wake)  are  likely  biased  by  probe- 
induced  oxidation  of  NO,  elevated  levels  of  NO2  In  areas  of  rapid  flame  quench 
(e.g.,  the  flame  front)  are  likely  real  and  not  artifacts  of  the  probe. 

Mainstream  and  Jet  Velocities.  The  major  effect  of  changing 
mainstream  and  jet  velocity  Is  to  change  the  size  of  the  recirculation  and 
wake  region.  This  Is  demonstrated  In  the  present  study  by  Independently 
Increasing  the  mainstream  velocity  (O^j)  and  decreasing  the  jet  velocity 
(Uj).  The  effect  of  either  Is  to  decrease  the  size  of  the  recirculation  and 
wake  regions. 

The  visual  appearance  of  the  flame  for  the  base  case  and  two 
variations  on  the  base  case  Is  shown  in  Figure  15c.  Both  the  penetration  of 
the  jet  and  the  radial  propagation  in  the  wake  are  restricted  by  Increasing 
the  mainstream  velocity  or  by  decreasing  the  jet  velocity.  This  is  confirmed 
by  the  detailed  temperature  maps  presented  In  Figure  15d. 

A  decrease  In  the  size  of  the  recirculation  and  wake  reaction  zones 
produce  a  net  reduction  In  the  residence  time  and,  hence,  a  net  reduction  In 
the  NOjj  production  (Figure  15e). 

Mainstream  and  Jet  Equivalence  Ratios.  The  mainstream  equivalence 
ratio  Is  the  dominant  variable  controlling  the  heat  release  and, 

ultimately,  the  pollutant  emission.  The  effect  on  heat  release  Is  shown  In 


Figure  15f. 


HC,  CO,  T.  As  Indicated  by  HC  exit  plane  profiles  (Figure  15f),  the 

wake  reaction  propagates  further  radially  as  the  mainstream  mixture  Is 

enriched.  This  Is  a  consequence  of  the  decrease  In  air  dilution  as  A  Is 

n 

Increased  from  0.6  to  1.0.  The  Increase  In  wake  reaction  as  (b  Is  enriched 

^m 
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from  1.0  to  1.2  Is  attributed  to  an  Increased  availability  of  hydrocarbon 

radicals.  Peak  flame  velocities  for  propane-air  flames  generally  occur  at 

equivalence  ratios  rich  of  stoichiometric  (Fristrom  and  Westenberg,  1965). 

The  highest  temperatures  occur  for  the  base  case  (A  -  1.0).  Peak 

m 

temperatures  are  about  300®K  lower  for  -  1.2  and  500*K  lower  for  0.8. 

Carbon  monoxide  concentrations  Increase  with  mainstream  equivalence 

ratio  as  the  amount  of  available  oxygen  to  oxidize  CO  to  CO2  decreases.  For 

all  conditions,  the  temperature  is  sufficient  for  the  oxidation  to  occur.  For 

equivalence  ratios  of  0.8  to  1.0,  peak  concentrations  correspond  to  the 

location  of  the  wake  reaction  front.  Inside  the  front,  CO  is  oxidized  to 

CO2.  This  accounts  for  the  Increase  in  temperature  In  the  wake.  Ahead  of  the 

front,  CO  diffuses  Into  the  cold  reactant  gases. 

At  «  1.2  ,  the  absence  of  oxygen  In  the  wake  results  In 
n 

relatively  constant  CO  concentrations  and  an  absence  of  a  distinct  CO  peak  at 
the  flame  front. 

For  all  cases  except  <{1^  ••  0.6,  temperatures  are  fairly  constant 

within  the  wake  reaction  zone  and  drop  at  the  radially  propagating  wake 

reaction  front.  For  ■■0.6,  the  temperatures  drop  immediately  adjacent  to 

m 

the  axis,  and  HC  concentrations  remain  elevated  while  CO  concentrations  fall. 

Instead  of  rise,  when  proceeding  from  the  axis  to  the  chamber  wall.  This 

suggests  that  reaction  In  the  wake  Is  suppressed  and  CO  formation,  for 

example.  Is  restricted  to  the  recirculation  region  upstream  with  radial 

diffusion  In  the  wake.  Note  that  CO  concentrations  for  A  -  0.6  are  not 

n 

appreciably  lower  than  for  ^  -  0.8.  In  fact,  near  the  jet  tube, 

m 

concentrations  are  lower  for  A  >■  0.8.  The  additional  oxygen  available 

Q 

at  ^  >  0.6  Is  offset  by  lower  temperatures, 
n 

Varying  the  jet  stream  equivalence  ratio  ((^  )  allows  determination 


of  the  effect  of  recirculation  zone  mixture  ratio.  The  effect  on  exit  plane 
profiles  of  HC  and  temperature  is  pronounced  only  at  (b  >■  0.6.  Higher  jet 

HI 

equivalence  ratios  result  in  lower  HC  concentrations  near  the  jet  tube  wall 
and  higher  temperatures,  and  this  effect  diminishes  as  distance  from  the  jet 
tube  wall  increases.  This  is  attributed  to  higher  temperatures  in  the 
recirculation  zone  that  result  from  recirculation  zone  mixture  ratios  closer 
to  stoichiometric.  This  effect  is  not  as  pronounced  for  the  other  equivalence 
ratios  because  each  effectively  sustains  a  fully  developed  reaction  in  the 
wake. 

The  effect  of  jet  equivalence  ratio  on  carbon  monoxide  is  noticeable 

only  near  the  jet  tube.  For  (J)^  ■  0.6,  richer  jet  mixtures  result  in  lower  CO 

emissions  because  of  higher  temperatures  in  combination  with  the  elevated 

concentration  of  oxygen.  This  same  trend  occurs,  but  to  a  much  lesser  extent, 

for  *  -  0.8.  At  (^  -  1.2  and  1.2  this  trend  is  reversed.  For  these  cases 
m  m 

oxygen,  and  not  temperature,  limits  the  CO  oxidation. 

NO  .  The  NO  profiles  are  presented  in  Figure  15g.  Peak 

X  X 

concentrations  are  highest  for  the  base  case  “  1.0).  The  concentrations 

are  slightly  lower  at  i  -1.2.  For  the  leaner  cases,  there  is  a  significant 

m 

drop  in  concentrations. 

For  all  cases  the  shape  of  the  NO^  exit  plane  profile  is  similar, 

decreasing  almost  linearly  from  the  jet  to  the  combustor  wall.  The  shape 

indicates  that  most  of  the  NO^  is  formed  in  the  recirculation  zone,  and 

diffuses  by  turbulent  transport  downstream.  This  is  confirmed  in  the  detailed 

flow  maps  presented  in  Figure  15g.  These  trends  correspond  well  with  the 

trends  observed  for  temperature  (Figure  15f)  reflecting  the  temperature 

dependence  of  NO^  formation  reactions. 

Note  that  the  ^  -  1.2  exit  plane  profiles  (Figure  15g)  intersect 
m 
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the  *1.0  profiles  with  higher  concentrations  near  the  chamber  wall. 

HI 

This  is  attributed  to  the  additional  production  of  in  the  larger  wake 
reaction  zone  associated  with  the  rich  mainstream.  Although  the  recirculation 
zone  Is  larger  as  well ,  the  data  Indicated  that  NO^  production  In  the 
recirculation  zone  Is  not  Increased,  a  consequence  of  suppressed  oxygen 
availability  and  temperature. 

Jet  equivalence  ratio  ((>^)  directly  Impacts  both  the  mixture  ratio 

and  temperature  of  the  recirculation  zone.  As  a  result,  the  effect  of  jet 

equivalence  ratio  on  NO^  production  is  predictable.  Production  of  NO^  is 

Increased  with  jet  enrichment  for  lean  mainstream  ■  0.6,  0.8),  decreased 

with  jet  enrichment  or  jet  leaning  for  stoichiometric  mainstream  “  1*0), 

and  Increased  with  a  lean  jet  or  decreased  with  a  rich  jet  for  a 

rich  (a  ••1.2)  mainstream, 
n 

The  effect  of  mainstream  equivalence  ratio  on  the  NO/NOj^  ratio  is 

shown  in  Figure  15h  for  a  stoichiometric  jet.  (Other  jet  equivalence  ratios 

are  omitted  for  clarity.)  The  rapid  drop  in  the  NO/NO^  ratio  occurs  at  the 

flame  front  for  each  of  the  cases  (6  >  1.2,  1.0,  0.8)  wherein  a  wake  reaction 

la 

was  supported.  The  low  N0/N0_  ratio  for  -  0.6  Is  attributed  to  the  quench 

X 

zone  surrounding  the  hot  recirculation  zone  In  the  absence  of  a  wake  reaction. 

Emission  Indexes.  The  emission  Indexes  for  NO^^,  CO,  and  HC  are 

presented  In  Figure  151  as  a  function  of  mainstream  equivalence  ratio  (A  ). 

m 

The  parameters  are  jet  equivalence  ratio  ((^j)  and  mainstream  velocity  (U^). 

The  procedure  used  to  compute  the  emission  Index  Involved  correcting  the  data 
for  water  vapor  In  the  exhaust,  calculating  the  area-averaged  exit  plane 
concentrations  and  the  area-averaged  mass  emission,  and  taking  the  ratio  of 
the  mass  emission  to  the  fuel  mass  Input. 

The  emission  Index  data  reflect  the  observations  derived  from  the 
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detailed  results  above,  and  place  the  performance  of  the  combustor  into  a 

practical  perspective.  For  example,  NO^  emissions  are  highest 

at  <b  ■  1.0,  and  are  only  slightly  lower  at  1.2.  Values  are  50  percent  lower 
m 

at  u  ■  0.8  than  at  stoichiometric,  and  those  at  0.6  are  from  5  to  15  times 
in 

lower  than  those  for  the  base  case  (depending  on  jet  equivalence  ratio).  The 
high  temperatures  that  favor  NO^  formation  also  favor  hydrocarbon  oxidation. 

In  the  absence  of  a  developed  wake  reaction  at  4)^  »  0.6,  60  percent  of  the 

fuel  Is  emitted  unburned.  As  suggested  by  the  analysis  of  the  detailed  flow 
maps,  jet  equivalence  ratio  has  a  minimum  Impact  on  HC  and  CO  emission,  but 
does  affect  NO^  emission. 

Summary .  Aerodynamic  flameholding  Is  shown  to  provide  control  over 
pollutant  formation  and  flame  stability  by  direct,  reactant  Injection  Into  the 
recirculation  zone.  The  data  presented  provide  Insight  into  the  performance 
of  the  reverse  jet  flameholder,  give  guidance  for  practical  application,  and 
establish  a  data  base  for  future  testing  of  elliptic  codes. 

The  reverse  jet  combustor,  chemically  and  aerodynamlcally,  consists 
of  two  distinct  regions:  the  recirculation  zone  and  the  wake.  The  mainstream 
and  jet  velocity  Influence  the  size  of  these  two  regions,  while  the  mainstream 
and  jet  mixture  ratios  affect  the  overall  chemistry  and  heat  release. 

The  Influence  of  the  jet  on  the  emission  of  NO^  Is  one  of  the  more 
Interesting  characteristics  of  the  flameholder.  Jet  changes  which  reduce  the 
size,  lower  the  temperature,  and  decrease  the  residence  time  are  favorable  to 
the  reduction  of  emission.  For  example,  an  enrichment  of  jet  equivalence 
ratio  affects  recirculation  zone  temperature  and  mixture  ratio,  and  will 
produce  either  an  Increase  or  decrease  in  the  net  emission  of  NO^  depending  on 
the  mainstream  equivalence  ratio.  Such  changes  do  not  significantly  affect  HC 
and  CO  emissions.  A  primary  benefit  from  an  enriched  jet  Is  to  extend  the 
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lean  blow-off  limit  and  thereby  maintain  combustor  stability  simultaneous  with 
a  reduction  in  the  emission  of  NO^.  However,  a  practical  limit  exists  beyond 
which  the  emission  of  unburnt  fuel  is  excessive.  In  the  present  experiment, 
this  limit  occurred  at  a  mainstream  equivalence  ratio  of  approximately  0.8. 

Detailed  information  is  available  in  a  M.S.  Thesis  (McDannel,  1979), 
and  a  Ph.D.  Dissertation  (Wuerer,  1978).  One  paper  has  been  published 
(McDannel,  Peterson,  and  Samuelsen,  1982). 

4.2.4  Dilute  Swirl  Combustor  (COLD  Experiments) 

Velocity  Data.  Axial  and  tangential,  mean  and  rms  velocity 
measurements  were  made  using  the  two-color  laser  anemometry  (LA)  system  shown 
in  Figure  3b.  Inlet  velocity  profiles  were  measured  as  closely  as  possible  to 
the  upstream  entry  plane.  Tangential  velocity  (w,  w^^g)  was  measured  1  mm  and 
axial  velocity  (u,  n^^)  was  measured  5  mm  downstream  of  the  entry  plane. 
Radial  velocities  were  not  measured  but  were  estimated  based  on  the 
extrapolation  of  streamlines  (calculated  from  measured  data)  to  the  entry 
plane.  These  inlet  velocity  profiles  are  presented  in  Figure  16a  for  the 
baseline  case  (15  m/s,  swirl/dilution  =  1.0,  COLD,  (|)  =  0.2).  The  spike  in 
the  inlet  tangential  velocity  profile  is  a  result  of  the  12  discrete  swirl 
channels.  Radial  velocity  measurements  across  the  full  combustor  diameter 
indicated  that  this  spike  was  a  symmetric  occurrence. 

Radial  and  centerline  profile  plots  of  mean  and  rms,  axial  (u,  n^^^) 
and  tangential  (w,  w^^^^g)  velocities  are  presented  in  Figure  16b  for  the 
baseline  case,  a  flow  condition  corresponding  to  a  15  m/s  reference  velocity 
and  equivalence  ratio  of  0.2.  Axial  values  appear  on  the  top  half  of  each 
radial  profile  plot  and  tangential  values  on  the  bottom.  Radial  profiles  were 
measured  at  seven  axial  stations  for  each  run  condition;  for  clarity  only  four 
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a)  Inlet  Profiles 


Figure  16.  DlluCe  Swirl  Combustor  COLD  Flowfleld  Data 
( 15  m/s,  COLD,  i  -  0.2) 


are  shovra  In  the  figures. 

Along  the  centerline,  tangential  velocities  go  to  zero  so  the  rms 
fluctuations  can  be  thought  of  as  rather  than  The  expansion  of  the 

swirling  air  Is  shown  to  form  a  strong  zone  of  backmlxlng  or  recirculation 
extending  about  10  cm  or  1.3  combustor  diameters  downstream.  The  peak  reverse 
flow  velocities  are  a  significant  fraction  (about  502)  of  the  bulk  velocity 
(15  m/s).  Proceeding  downstream,  the  radial  profiles  of  mean  axial  velocity 
show  the  extent  of  the  recirculation  zone  and  the  trend  towards  a  pipe-flow 
profile. 

The  radial  profiles  of  mean  tangential  velocity  show  that  the 
Initially  strong  concentrated  swirling  air  mixes  with  the  dilution  air  and 
forms  a  more  uniform  distribution  of  swirl  downstream  with  a  much  lower 
maximum  velocity.  It  can  be  seen  that  the  dilution  of  the  swirl  strength 
corresponds  with  closure  of  the  recirculation  zone.  At  the  upstream  end  of 
the  recirculation  zone  the  step  In  the  radial  profile  of  tangential  velocity 
Is  evidence  that  the  swirl  near  the  centerline  at  this  axial  station  Is 
convected  upstream  from  a  downstream  location.  At  the  most  downstream  radial 
profile,  the  axial  velocities  are  seen  to  rapidly  approach  a  "well  developed" 
turbulent  pipe  flow  profile  while  the  tangential  velocities  are  somewhat 
slower  In  approaching  a  "solid  body  rotation"  form. 

Radial  profiles  of  turbulence  Intensity  show  peaks  of  the  regions  of 
maximum  shear.  At  the  upstream  regions  of  the  flowfleld,  the  axial  rms 
velocity  (Uj^g)  Is  generally  larger  than  the  tangential  rms  velocity  (w^^g). 
However,  Isotropy  seems  to  be  a  reasonable  engineering  assumption  for  the  case 
since  w^_  Is  generally  within  202  of  u^_  except  In  the  regions  where  (1)  the 
non-swlrllng  fuel  Jet  Issues  (Illustrated  by  the  peak  In  the  u  velocity 
curve),  and  (2)  the  non-swlrllng  dilution  air  enters. 
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Using  the  continuity  equation,  stream  functions  t|i  can  be  obtained  by 
radially  Integrating  the  profiles  of  mean  axial  velocity: 

^rudr  where  r  «  radial  location  of  interest 

I 

rudr  R  >  outer  radius  of  test  section 

Integration  at  each  profile  allows  streamlines  to  be  drawn  through  points  of 
constant  i|i.  These  streamlines,  shown  In  Figure  16c  Illustrate  the  form  of 
the  "time-averaged”  flowfleld. 

Figure  16d  Is  a  photograph  of  streak  lines  created  by  laser 
Illuminated,  neutrally  buoyant,  helium-filled  bubbles  that  were  Injected  Into 
the  DSC  to  provide  flow  visualization.  Comparison  of  these  streak  lines  to 
the  time-averaged  streamlines  shows  agreement  In  the  length  and  extent  of  the 
recirculation  zone.  The  streak  lines  depict  a  high  degree  of  mixing  and 
disorder  in  the  vicinity  of  the  recirculation  zone. 

Modeling.  The  TEACH  finite  difference  code  with  the  k-e  turbulence 
model  (Gosman  and  Pun,  1974)  was  modified  to  Incorporate  swirl  following  Lilly 
and  Rhode  (1982)  and  used  to  predict  the  DSC  flowfleld  characteristics.  The 
object  was  to  demonstrate  the  Isothermal  (COLD)  performance  and  ease  of 
applicability  of  this  commonly  used  turbulence  model  on  a  swirling  complex 
flow. 

The  code  was  applied  to  the  boundary  conditions  of  the  DSC  flowfleld 
for  the  baseline  case  (15  m/s,  ■  0.2,  COLD)  without  difficulty.  Calculated 
variables  were  the  three  components  of  velocity  (u,v,w),  pressure  (p), 
turbulent  kinetic  energy  (k),  and  dissipation  rate  (e).  The  code  was  run 
with  the  standard  k-e  turbulence  model  and  radial  profiles  of  the  predicted 
mean  axial  and  tangential  velocities  are  compared  with  measured  results  In 
Figure  16e.  The  code  predicts  a  much  slower  radial  diffusion  of  both  u  and  w 
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velocity  gradients  than  Is  apparent  from  the  experimental  results.  It  falls 
to  predict  the  strength  and  size  of  the  recirculation  zone.  Similar  results 
have  been  reported  by  others  when  modeling  swirling  flows  and  the  differences 
are  often  attributed  to  streamline  curvature  (e.g.,  Srlnivasan  and  Hongla, 
1980).  The  DSC,  as  a  result.  Is  an  attractive  configuration  In  which  to  test 
refinements  to  elliptic  codes  designed  to  Improve  performance  In  swirl 
dominated  flows  as  well  as  reacting  flows. 

Summary .  A  summary  of  this  work  is  Included  as  part  of  the  summary 
presented  at  the  end  of  the  next  section. 

4.2.5  Dilute  Swirl  Combustor  (HOT  Experiment) 

Velocity  Data.  Detailed  laser  anemometry  measurements  were  made  In 
reacting  flow  at  two  overall  equivalence  ratios  (0.2  and  0.1).  Centerline 
profiles  of  mean  and  rms  axial  (u,  u^^g)  velocities  appear  in  Figure  17a.  The 
recirculation  zone  for  the  ■  0. 1  case  appears  very  similar  to  that  for  the 
cold  flow  case  (Figure  16b)  with  three  notable  exceptions.  The  recirculation 
zone  In  the  case  of  reaction  Is  (1)  more  Intense  (higher  negative  velocities), 
(2)  more  compact  (shorter),  and  (3)  radially  wider.  Root-mean-square  velocity 
levels  are  approximately  50%  higher  for  the  reacting  case  except  In  the  center 
of  the  recirculation  zone  where  the  tms  levels  are  about  equal.  From  a  visual 
perspective,  the  flame  was  detached  from  the  face  of  the  fuel  delivery  tube 
and  was  uniformly  blue  with  no  Indication  of  soot  formation. 

As  the  equivalence  ratio  was  Increased  to  ^  -  0.2,  a  transition 
occurred  In  the  form  of  the  recirculation  zone.  Visually,  the  flame  became 
attached  to  the  upstream  fuel  tube  face  and  yellow  streaks  of  soot  were 
observed  emlnating  from  that  point  downstream  along  the  centerline  of  the 
combustor.  The  transition  Is  clearly  evidenced  by  the  centerline  profile  of 
u-veloclty  which  Indicates  the  presence  of  positive  velocities  just  downstream 
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and  Radial  Profiles  of  Axial  and  Azimuthal 


(15  m/s,  <^  -  0.1  and  0.2) 
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Figure  17.  Dilute  Swirl  Combustor  HOT  Flowfleld  Data  (Velocity) 
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d)  Radial  Profiles  of  the  Correlation  Coefficient  (15  m/st  <^  ■  0,2) 
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Figure  17.  (concluded) 


of  the  fuel  delivery  tube. 

The  effect  of  equivalence  ratio  is  evident  as  well  in  the  radial 
plots  of  axial  and  tangential  velocities  but  is  best  Illustrated  by  the 
streamlines  which  are  plotted  in  Figure  17b.  As  the  equivalence  ratio  is 
increased,  the  recirculation  zone  becomes  shorter  and  wider.  Eventually,  the 
recirculation  zone  becomes  so  wide  chat  not  all  the  conical  fuel  jet  is 
retained  in  the  outer  stream  layer.  Instead,  some  of  the  fuel  is  directed 
Inward  towards  Che  centerline  where  it  (1)  opposes  Che  upstream  momentum  of 
the  recirculating  flow  and  (2)  ignites  in  Che  relatively  oxygen  deficient  core 
of  Che  recirculation  zone  resulting  in  visible  soot  production. 

The  change  in  width  and  length  of  Che  recirculation  zone  associated 
with  reaction  is  a  result  of  gases  expanding  in  the  shear  layer  at  the 
boundary  of  Che  recirculation  zone  to  create  a  source  of  pressure  which  tends 
to  force  (1)  early  closure  of  the  recirculation  zone,  and  (2)  swirling  air 
Cowards  the  centerline.  As  Che  swirling  air  is  forced  towards  the  centerline, 
the  tangential  velocity  increases  (conservation  of  angular  momentum).  This 
increase  is  evident  by  comparing  the  tangential  velocity  (w)  profiles  of 
Figures  16b  and  17a.  When  this  Increased  tangential  velocity  is  convected 
upstream  through  the  recirculation  zone,  it  results  in  a  higher  radial 
pressure  gradient  that  tends  to  widen  the  zone. 

Turbulence  Intensity  levels  in  Che  vicinity  of  Che  recirculation 
zone  are  noticeably  higher  for  Che  ^  •  0,2  case  (in  comparison  to 
the  -  0.1  case),  whereas  they  are  about  equal  downstream.  Part  of  this 
Increased  level  is  attributed  Co  intermittency  in  Che  division  of  Che  fuel  Jet 
between  the  inside  and  outside  of  the  recirculating  flowfield. 

A  lower  reference  velocity  is  shown  in  Figure  17c  where  centerline 
profiles  of  and  velocities  for  the  7.5  m/s,  (^*0.2  reacting  case 


are  presented.  A  peak  in  both  the  and  velocities  occurs  at  the 

stagnation  point  of  the  recirculation  zone.  This  characteristic  has  been 
observed  by  others  (e.g.,  Owen  1976)  and  is  attributed  to  "form  intermittency” 
of  the  recirculation  zone.  Such  large  scale  intermittency  or  periodic 
behavior  of  complex  flows  has  been  documented  utilizing  flow  visualization 
techniques  such  as  high  speed  cinematography  (Bnim,  Ikioka,  and  Samuelsen 
1982)  and  time  resolved  emissions  techniques  (Roquemore,  Britton,  and  Sandhu, 
1982). 

It  is  also  observed  that,  at  downstream  locations  along  the 
centerline,  is  greater  than  u^g.  Examination  of  radial  profiles  of 

reveals  that  these  unusually  high  values  occur  only  near  the  centerline. 

These  excessive  levels  are  most  likely  produced  by  a  spiralling  action  of  the 
vortex  center  about  the  stationary  LA  probe.  This  spiralling  action  (a  type 
of  "form  intermittency"  or  "periodicity")  is  further  evidence  by  a  flattened 
probability  distribution  function  (pdf)  which  one  would  expect  with  a 
spiralling  center  versus  a  Gaussian  distribution  more  typical  of  "oure" 
turbulence.  Similar  high  turbulence  production  rates  have  been  observed  by 
Cheng  and  (1981)  when  performing  LA  measurement  on  a  premixed  turbulent 
flame.  They  attributed  this  to  periodic  movement  of  the  flame  sheet  through 
the  measurement  volume.  Hence,  the  "turbulence  intensity"  as  measured  by  a 
stationary  LA  probe  is  comprised  of  two  components:  (1)  rms  fluctuations 
intrinsic  to  fluid  mechanic  turbulence  production  and  (2)  form  intermittency 
of  the  complex  flowfleld  that  does  not  give  rise  to  significant  turbulent 
stress  production. 

Radial  plots  of  the  u'w'  correlation  coefficient  (C)  appear  in 
Figure  17d  (top  and  bottom)  for  the  baseline  (COLD)  case  along  with  the  axial 
(u,  top)  and  tangential  (w,  bottom)  velocity  profiles.  These  data  represent 


some  of  the  first  complex  swirling  flow  velocity  correlations  that  have  been 
obtained  by  direct,  simultaneous  measurement.  For  well  developed  boundary 
layer  flows,  the  coefficient  (C)  has  been  classically  measured  to  be  on  the 
order  of  0.4  (Llepmann,  1979).  In  complex  flows,  of  course,  the  direction  of 
the  coefficient  changes  as  does  the  direction  of  local  gradients.  For 
example,  where  both  the  u  and  w  gradients  with  respect  to  y  are  positive: 


then  the  correlation  (C)  is  positive;  when  both  are  negative,  the  correlation 
(C)  Is  positive.  However,  when  the  u  gradient  Is  positive  and  the  w  gradient 
is  negative,  or  vice  versa,  the  correlation  is  negative.  In  areas  where  there 
is  little  or  no  gradient  in  either  u  or  w,  the  correlation  goes  to  zero. 

The  peak  value  of  the  correlation  coefficient  measured  was 
approximately  0.25,  which  is  less  than  that  measured  In  boundary  layers 
(0.4).  However,  there  is  no  reason  to  believe  that  they  should  be  the  same 
since,  in  boundary  layers,  there  is  a  gradient  in  only  one  velocity 
component  (~  *  0)  and  not  both  as  in  the  present  case, 

3y 

(l.e.,  ^  0).  Smith  and  Giel  (1979)  obtained  peak  values  of 

the  u'w’  correlation  coefficient  of  0.23  via  a  similar  direct  method  of 
measurement  in  a  complex,  but  non-swirling,  flow  (dump  combustor). 

Indirect  measurements  of  u*w'  have  been  made  in  swirling  complex 
flow  geometries  under  Isothermal  conditions.  For  example,  Janjua,  et  al. 
(1982)  used  a  six  orientation  hot  wire  probe  in  an  expanding  swirl  combustor 
and  derived  the  various  flow  quantities  algebraically  based  on  time-averaged 
quantities  obtained  at  each  orientation.  The  peak  u*w'  correlation 
coefficient  calculated  from  their  results  was  on  the  order  of  0.38.  Fujll, 
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Eguchl,  and  Goml  (1981)  used  measurement  of  a  single  component  LA  system  made 
In  a  swirling  jet  at  various  angles  to  algebraically  determine  u'w'  .  The 
peak  coefficient  calculated  from  their  results  was  about  0.40.  Part  of  the 
difference  between  the  direct  and  Indirect  measurement  of  the  correlation 
coefficients  may  be  attributed  to  the  difference  in  flow  geometries.  For 
example,  the  level  of  "form  Intermittency"  In  the  present  complex  flow 
geometry  may  be  different  than  that  In  the  geometries  used  by  Fujll  and 
Janjua.  Since  the  stationary  LA  probe  Is  not  capable  of  distinguishing 
between  fluctuations  arising  from  (1)  true  turbulence  and  (2)  form 
Intermittency,  and  since  these  two  fluctuating  components  have  different  time 
scales  and  are  thus  not  well  correlated  to  one  another,  the  presence  of  the 
latter  (form  Intermittency)  effectively  dilutes  the  correlation  produced  by 
the  first  (true  turbulence). 

Reynolds  stress  plots  for  the  15  m/s,  (ji  •  0.2  reacting  flow  case 
appear  In  Figure  17d;  once  again  the  peak  correlation  coefficient  is  on  the 
order  of  0.25. 

Velocity/Temperature  Data.  As  a  next  step,  a  fine-wire  thermocouple 
probe  was  employed  to  obtain  temperature  measurements  simultaneously  with 
velocity. 

First,  the  results  of  a  probe  perturbation  study  are  presented  where 
the  peturbatlon  caused  by  the  use  of  probes  of  various  shapes  and  orientations 
Is  assessed.  Second,  the  effect  of  time  constant  variation  Is  considered 
followed  by  a  presentation  of  the  temperature  time  series  and  the  probability 
density  function  (PDF)  of  the  temperature.  The  section  concludes  with  the 
statistical  properties  of  the  combined  velocity  and  temperature  field.  A 
tabulated  data  base.  Including  measured  Inlet  conditions.  Is  available 
(Seller,  1983b). 


The  probe  perturbation  effects  are  both  small  scale  (local)  and 


large  scale  (global).  The  former  Is  associated  with  the  geometrical 
characteristics  of  the  probe,  especially  near  the  probe  tip,  while  the  latter 
Is  related  to  the  overall  shape,  size,  and  orientation  of  the  probe.  The 
local  effects  can  be  minimized  by  making  the  probe  tip  as  small  as  possible 
while  the  global  effects  can  be  minimized  by  probe  geometry  (including  overall 
size)  and  probe  orientation. 

Probe  perturbation  was  assessed  by  separately  placing  each  of  the 
probe  configurations  (Figure  4a)  In  the  flow  and  determining  statistical 
properties  of  the  velocity  field  in  the  presence  and  absence  of  the  probe.  In 
the  case  of  the  curved  probe,  different  orientations  of  the  probe  were  also 
assessed.  Representative  results  based  on  500  data  samples  are  presented  in 
Figure  18a.  The  probe  with  minimum  overall  perturbation  effects  is  identified 
by  the  solid  symbol.  In  the  central  region  of  the  flow,  the  curved  probe  has 
a  lower  perturbation  effect  than  that  of  the  straight  probe.  This  may  be  due 
to  the  occasional  occurrence  of  periods  of  reverse  flow  and/or  suppression  of 
the  form  intermlttency  observed  by  Brum  and  Samuelsen  (1982b). 

Based  on  this  perturbation  study,  probe  configurations  and 
orientations  were  selected  that  minimize  the  perturbation.  The  resultant 
uncertainty  due  to  probe  perturbation  were  less  than  25Z  throughout  much  of 
the  flow  field. 

The  effect  of  time  constant  variation  on  statistical  features  of  the 
flow  involving  the  temperature  signal  is  shown  in  Figure  18b.  The  mean 
temperature  is  not  significantly  affected  by  deliberately  over-  or  under¬ 
compensating  the  sensor  relative  to  the  value  obtained  with  the  sensor 
compensated  for  the  mean  flow  conditions.  For  example,  the  measured  value  of 
the  mean  temperature  varies  by  less  than  10%  for  a  ±  50%  variation  in  the  time 
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constant.  The  root-mean-square  temperature  Is  also  relatively  Insensitive  to 
variation  In  the  time  constant.  A  lOZ  variation,  has  less  than  a  5Z  effect  on 
the  measured  root -mean-square  temperature  which  suggests  that  over-  and  under¬ 
compensation  affects  that  portion  of  the  temperature  power  spectrum  that 
contributes  a  relatively  small  amount  to  the  root-mean-square  temperature.  In 
contrast,  a  20Z  variation  In  the  measured  axial  heat  flux  and  a  50Z  variation 
In  the  tangential  heat  flux  occurs  for  the  same  10%  under-compensation  In  the 
time  constant.  Over-compensation  has  relatively  little  effect  on  the  value  of 
the  heat  fluxes.  This  observation  Is  consistent  with  that  of  Yanagl  and 
Mlmura  (1981).  It  Is  clear  that  both  phase  shift  and  variations  In  the 
amplitude  response  that  are  correlated  with  the  temperature  and  velocity 
affect  the  value  of  the  heat  fluxes,  but  the  difference  between  the  effect  of 
over-  and  under-compensation  Is  difficult  to  explain. 

Representative  samples  of  the  compensated  temperature  signal  are 
presented  In  Figure  Idc  at  different  radial  positions  for  the  14  cm  axial 
station.  A  peak-to-peak  temperature  of  about  SOO’C  Is  found  at  the  centerline 
distributed  about  a  mean  temperature  of  1560°C.  The  fluid  on  the  centerline 
Is  well  mixed  with  no  evidence  of  unheated  dilution  air  reaching  the 
centerline. 

The  temperature  time  series  at  the  mld-radlus  position  Indicates 
that  the  temperature  signal  has  positive  skewness  due  to  the  presence  of 
relative  short  duration,  high  temperature  peaks.  The  average  rate  of 
occurrence  of  hot  particles  (defined  as  the  average  rate  at  which  the 
temperature  signal  crosses  a  threshold  temperature)  Is  42  Hz  for  a  threshold 
of  lOOO'C.  The  crossings  are  not  periodic  but  appear  to  be  distributed 
randomly  In  time,  and  are  evidence  that  a  periodic  structure  Is  not  associated 
with  the  passage  of  the  high  temperature  fluid  particles.  Power  spectral 


neasurements  of  the  temperature  signal  show  a  relative  Increase  in  the 
“energy"  of  the  temperature  signal  at  a  frequency  of  90  Hz.  This  corresponds 
to  the  frequency  of  longitudinal  oscillation  of  the  recirculation  zone 
observed  In  high  speed  photography  (Brum  and  Samuelsen,  1982b). 

Hear  the  wall,  the  temperature  signal  Is  dominated  by  high 
amplitude,  short  duration,  temperature  spikes  which  correspond  to  a  positive 
skewness.  The  average  rate  of  occurrence  of  hot  particles  crossing  a 
threshold  of  400”C  Is  25  Hz.  Again,  the  crossings  are  not  periodic.  The 
minimum  temperature  corresponds  to  about  40" C.  The  data  obtained  at  the  24  cm 
station  are  similar.  Indicating  that  the  sheath  of  dilution  air  near  the  wall 
of  the  combustor  remains  Intact  through  the  length  of  the  combustor. 

Representative  probability  density  functions  (PDF)  are  presented  In 
Figure  18d  for  two  axial  locations,  14.0  and  24.0  cm.  The  data  are  normalized 
by  the  maximum  adiabatic  temperature,  T^j  -  1996"C.  The  PDF's  show  positive 
skewness  at  locations  well  displaced  from  the  centerline  (r/R  >  0.6).  This 
was  obtained  at  all  four  axial  stations  and  Indicates  that  cool,  unreacted  gas 
predominates  In  this  region  of  the  flow  throughout  the  length  of  the 
combustor.  In  addition,  the  PDF's  are  narrow  In  this  outer  region  of  the  flow 
which  corresponds  to  low  rms  temperatures  and  Indicates  that  only  a  limited 
amount  of  fluid  from  the  hot  recirculation  zone  mixes  with  the  cool  dilution 
air. 

At  the  14.0  cm  station,  the  PDF's  are  very  broad  at  the  intermediate 
radial  locations  of  r/R  -  0.4  and  r/R  •  0.3  with  peak-to~peak  temperatures  on 
the  order  of  1100"C,  Indicating  the  alternating  presence  of  cool  and  hot  fluid 
parcels.  Maximum  rms  temperatures  are  also  found  at  these  radial  locations. 

At  r/R  <0.2,  the  gas  Is  mainly  high  temperature  (1160"C  to  1750”C)  combustion 
products  with  no  evidence  of  unheated  dilution  air.  At  the  24  cm  station,  the 


data  provide  evidence  of  additional  mixing  between  the  hot  products  of 
combustion  and  cool  dilution  and  unreacted  swirl  air.  The  PDF's  at  r/R  <  0.2 
are  broadened  relative  to  the  upstream  14  cm  station,  and  the  peak 
temperatures  of  the  PDF's  at  the  Intermediate  radial  locations  (r/R  •  0.3, 

0.4)  are  reduced. 

Velocity  data  obtained  with  the  thermocouple  probe  In  place  (Figure 
18e)  Indicate  that  the  statistical  features  of  the  flow  field  remain 
essentially  unchanged  from  those  found  In  the  absence  of  a  probe  (cf..  Figure 
17a). 

The  mean  and  rms  temperature,  the  axial  velocity,  and  the  normalized 
axial  heat  flux  are  shown  In  Figure  18f.  The  mean  temperature,  at  all  axial 
locations,  has  a  maximum  on  the  centerline.  The  maximum  mean  temperature  In 
the  recirculation  region  at  x  >  2.0  and  7.0  cm  Is  nearly  constant  and  equal  to 
nsO’C  which  Is  approximately  200*C  less  than  the  maximum  adiabatic  flame 
temperature  (1996‘C).  Radiation  losses  are  estimated  to  account  for 
approximately  90*’C  of  the  200*’C.  The  nearly  constant  mean  temperature  In  the 
recirculation  region  Indicates  that  the  fluid  In  this  region  Is  relatively 
well  mixed.  The  peak  Instantaneous  temperature  recorded  was  as  high  as 
2094 "C.  The  Infrequent  occurrence  of  measured  temperature  greater  than  the 
maximum  adiabatic  temperature  occurred  less  than  \X  of  time  and  Is  due  to 
slight  over-compensation  of  the  signal.  At  x  «  2.0  and  7.0  cm  for  r/R  >  0.7, 
the  mean  temperature  Is  less  than  ISO'C  and  consequently  the  flow  corresponds 
to  slightly  heated  dilution  air. 

Downstream  of  the  recirculation  zone  (at  x  -  14.0  and  24.0  cm),  the 
mean  temperature  profile  exhibits  a  maximum  of  1400”C  on  the  centerline  with  a 
monotonlc  decay  to  a  minimum  near  the  wall  of  125 **0. 

The  axial  heat  fluxes  at  the  14.0  and  24.0  cm  stations  are  negative 


In  the  central  region  of  the  flow  (r/R  <  0.3)  and  are  either  relatively  small 
or  slightly  positive  In  the  outer  region  of  the  flow  (r/R  >  0.4).  The 

negative  sign  of  the  axial  heat  flux  In  the  central  region  of  the  flow 
Indicates  the  presence  of  a  relatively  large  number  of  fluid  particles  with 
high  temperature  (relative  to  the  mean  temperature)  and  low  velocity  (relative 
to  the  mean  velocity).  These  fluid  particles  probably  originate  In  the 
recirculation  zone.  The  growth  of  the  region  of  negative  axial  heat  flux  In 
the  downstream  direction  Is  associated  with  the  mean  tangential  velocity  which 
preferentially  transports  low  temperature »  high  density  fluid  away  from  the 
central  region  and  high  temperature,  low  density  fluid  toward  the  centerline. 

The  mean  tangential  velocity  at  the  2.0  cm  station  (Figure  18g)  Is 
relatively  low  In  the  recirculation  region,  with  the  peak  mean  tangential 
velocity,  as  expected.  In  line  with  the  outer  circumference  of  the  swirl 
vanes.  This  sharp  peak  Is  quickly  suppressed  (x  «  7.0  cm)  by  the  Interaction 
with  the  non-swlrl  dilution  air.  The  tangential  heat  flux  Is  very  small  at 
the  X  >■  2.0  cm  station.  Downstream  (x  >  14.0  and  24.0  cm),  the  values  of  the 
tangential  heat  flux  Increase  and  exhibit  relatively  sharp  peaks  at  radial 
positions  (r/R  -  0.2  and  0.4  respectively)  corresponding  to  the  maximum  values 
of  the  mean  tangential  velocity. 

The  mean  temperature  In  the  central  region  at  x  -  2.0  and  7.0  cm  Is 
so  high  that  the  overheating  of  the  probe  required  In  order  to  determine  the 
time  constant  caused  the  probe  to  fall.  Therefore,  rms  temperature  and  axial 
and  radial  heat  fluxes  could  not  be  determined  In  this  region. 

Conditioned  sampling  techniques  were  applied  to  the  data  In  which  a 
temperature  range  Is  specified  and  the  statistical  quantities  are  obtained  for 
those  data  points  falling  within  the  specified  range.  The  results  verify  that 
the  central  region  of  the  combustor  downstream  of  the  recirculation  zone  Is 


dominated  by  parcels  of  high  temperature,  low  velocity  fluid.  At  the  14  cm 
station  near  the  centerline,  for  example,  less  than  one  percent  of  the  fluid 
parcels  have  simultaneously  high  temperature  and  high  velocity. 

Summary.  Non-intruslve  two-color  laser  anemometry  has  been  used  to 
perform  a  flowfleld  characterization  of  a  model  complex-flow  combustor.  This 
Is  an  Important  first  step  toward  the  establishment  of  a  "pool"  of  benchmark 
data  on  a  common  complex  flow  test  bed  geometry  that  can  be  used  to  (1) 
evaluate  turbulence  models,  and  (2)  enhance  the  general  understanding  of 
turbulent  processes  In  practical  complex  flows. 

A  technique  was  used  to  obtain  a  direct  measurement  of  Reynolds 
stress  (u'w' ) ;  peak  correlation  coefficients  were  on  the  order  of  0.25  which 
is  less  than  that  reported  by  others  (0.40)  using  indirect  measurement  methods 
on  somewhat  different  complex  flow  geometries. 

Root-mean-square  velocity  results  indicate  that  for  this  flow  u-w 
Isotropy  is  a  reasonable  engineering  assumption.  Except  In  a  few  regions 
(e.g.,  the  area  of  the  fuel  Jet),  the  values  of  and  w^^  are  generally 

within  20%. 

The  effect  of  reaction  was  to  Increase  turbulence  Intensity  levels 
and  the  equivalence  ratio  was  observed  to  have  a  significant  effect  on  the 
size  and  form  of  the  recirculation  zone.  Fuel  Injection  techniques, 
investigated  In  previous  work  (Brum  and  Samuelsen,  1982a),  Indicate  that 
"tailoring"  fuel  distribution  to  the  recirculation  zone  shape  Is  Important  In 
controlling  soot  production.  In  the  present  work,  quantitative  flowfleld 
measurements  verify  that  the  Increased  visible  soot  production  at  an 
equivalence  ratio  of  ■  0.2  corresponds  with  a  deflection  of  part  of  the 
fuel  jet  towards  the  relatively  oxygen  deficient  core  of  the  recirculation 
zone.  Fuel  injection  techniques  which  Insure  the  placement  of  fuel  into  the 
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outer  shear  layer  of  the  recirculation  zone  reduce  the  dependence  of  sooting 
propensity  on  equivalence  ratio. 

There  is  evidence  of  "form  intermit tency”  of  the  recirculation  zone 
which  can  take  one  of  two  forms:  (1)  a  spiralling  action  of  the  vortex 
center,  and  (2)  a  fluctuating  interaction  between  the  fuel  Jet  and  the 
recirculation  zone.  The  magnitude  of  this  intermit tency,  particularly  the 
latter,  is  higher  reacting  flew  and  a  function  of  equivalence  ratio.  When 
measured  with  a  stationary  LA  probe,  this  form  Intermlttency  manifests  Itself 
in  the  form  of  increased  turbulence  intensity  levels;  this  component  tends  to 
"dilute"  the  Reynolds  stress  correlation  coefficient. 

Compensated  fine-wire  thermocouples  can  be  used  to  measure  time- 
resolved  temperature  along  with  axial  and  tangential  velocity  in  a  complex 
flow,  swirl-stabilized  combustor  but  precautions  are  required,  especially  with 
respect  to  probe  perturbation  and  survivability.  Different  probe 
configurations  and  different  probe  orientations  must  be  used  in  the  various 
regions  of  the  flow  to  minimize  probe  perturbation.  The  sensor  size  required 
to  survive  in  the  flow  depends  on  the  mean  temperature  and  velocity  intensity 
and  must  be  Increased  in  high  temperature,  highly  turbulent  regions.  In  the 
high  temperature  core  of  the  recirculation  zone,  the  time  constant  cannot  be 
accurately  measured  due  to  probe  failure  during  the  overheating  process. 

The  mean  temperature  is  nearly  constant  in  the  recirculation  zone 
with  peaks  approximating  the  maximum  adiabatic  flame  temperature.  Further 
downstream,  mixing  of  hot  products  with  unreacted  swirl  and  dilution  air  takes 
place  and  a  gradient  in  tiie  mean  temperature  develops  in  the  central  region. 
Axial  heat  fluxes  in  the  central  region  are  negative,  reflecting  a  large 
population  of  high  temperature,  low  velocity  fluid  from  the  recirculation 
region.  The  swirl  velocity  is  shown  to  increase  this  population  by  inducing  a 


transport  of  low  density,  high  temperatures  fluid  to  the  axis. 

Detailed  information  is  available  on  the  two-component  Laser 
Anemometry  data  in  a  Fh.D.  Dissertation  (Brum,  1983).  One  paper  has  been 
presented  (Samuelsen  and  Brum,  1982b)  and  a  second  has  been  invited  for 
presentation  at  the  Symposium  on  Experimental  Measurements  and  Techniques  in 
Turbulent  Reactive  and  Non-Reactlve  Flows,  ASME  105th  Winter  Meeting  (Brum  and 
Samuelsen,  1984b).  A  tabulated  data  set  is  available  (Brum,  1983). 

Four  papers  have  been  presented  on  the  combined  velocity-temperature 
measurements  (Brum,  et  al.,  1983;  Seller,  et  al.,  1983;  Samuelsen,  LaRue,  and 
Seiler,  1984;  LaRue,  Samuelsen,  and  Seiler,  1984)  and  one  paper  has  been 
accepted  for  publication  (LaRue,  Samuelsen,  and  Seller,  1984).  Detailed 
information  is  available  in  a  M.S.  Thesis  (Seller,  1983a).  A  tabulated  data 
base  is  provided  in  Seiler  (1983b). 

4.3  Element  C;  Supplemental  Studies 
4.3.1  LA  Seeder  Development 

Introduction.  In  non-premixed  combustion  systems,  fuel  and  oxidant 
a'  Injected  into  the  reaction  chamber  via  two  or  more  separate  stream. 
Velocity  measurements  by  laser  anemometry  within  the  chamber  may  be  biased  by 
individual  streams  unless  the  seeding  concentrations  of  both  streams  are 
uniform  and  constant.  The  present  study  addressed  the  development  of  a  seed 
generation  technique  for  combusting  flows  and  the  evaluation  of  its 
performance  in  providing  both  a  steady  rate  of  injection  and  a  uniform 
volumetric  seed  concentration  for  a  combustion  chamber  featuring  multiple 
stream  injection. 

Particle  Generator.  The  particle  generation  seeding  system  is 
schematically  depicted  in  Figure  19a.  Air  is  delivered  at  line  pressure 
(~60  pslg)  to  a  dryer  (Wilkerson  Model  X03-02-000)  for  the  removal  of  residual 
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a)  Seeding  System 


b)  Particle  Rate  Measurement  System 


Figure  19.  LA  Seeder  Results 


moisture.  The  air  Is  then  channeled  through  a  high-quality  filter  (Balston  DQ 
Grade),  capable  of  90Z  retention  of  0.6  von  >  remove  undesirable  residue. 

The  air  pressure  Is  regulated  (Norgren  Type  R-11),  usually  20-30  pslg,  to 
provide  atomization  conditions  specified  by  the  atomizer  (RETEC  X70).  The 
atomization  air  and  the  secondary  air  are  metered  by  rotameters  (DWYER  Hodels 
RICB-51  and  RMB-54),  ordinarily  In  10:1  secondary  to  atomization  air 
proportions.  A  suspension  of  metal  oxide  (aluminum  oxide,  Al202>  In  the 
present  case)  and  deionized  water  Is  agitated  by  a  magnetic  stirrer  to  provide 
a  uniformly  mixed  suspension.  The  suspension  Is  mixed  In  proportions,  usually 
4  g/1,  that  minimize  the  probability  of  having  multiple  particles  per  drop. 

The  atomizer  produces  an  aerosol  with  a  mass  median  diameter  (MMD)  of 
5.1  ym  with  a  geometric  standard  deviation  of  2  ym  .  After  atomization,  the 
aerosol  Is  carried  to  a  mixing  chamber.  The  secondary  air  Is  Injected 
tangentially  Into  the  mixing  chamber  to  Induce  mixing  and  decrease  the 
humidity  of  the  seeded  flow.  The  capability  of  charge  neutralizing  the 
aerosol  Is  also  provided  by  passing  the  aerosol  through  a  radiation  flux 
produced  by  a  10-mC  vial  of  Krytpon-85  gas  (3M  Model  3B4G). 

System  Evaluation.  Light  scattering  refractory  particles  of 
sufficient  size  to  provide  adequate  light  but  small  enough  to  follow  the  flow 
are  required  for  laser  measurements.  Thus,  1.0  ym-dlam  ^2^3  P^^^^^cles  were 
selected  for  the  present  case  to  provide  adequate  scattering  and 
representation  of  fluid  velocity. 

Al20^  particles  nominally  sized  at  1.0  ym  and  1.5  ijn  were  used.  The 
actual  size  of  generated  particles  was  determined  by  a  scanning  electron 
microscope  (SEM).  Micrographs  established  that  the  nominal  1.0  ym  amorphous 
Al20^  particles  were  substantially  smaller,  approximately  0.1  •  The 

significant  variation  from  the  quoted  size  In  the  l.O-im  ^2^3 


to  the  morphology  o£  the  Al20^.  Amorphous  AI2O2  is  composed  of  agglomerations 
of  ~0.1-ym  particles  having  an  average  agglomerated  diameter  of  1  ym.  When 
Al20^-H20  suspensions  are  atomized  and  the  moisture  evaporated,  the  resultant 
particles  are  considerably  smaller  than  the  quoted  size.  In  hot  flow  cases, 
the  sudden  heating  breaks  apart  the  remaining  agglomerations  that  are  of 
sufficient  size  to  be  useful,  causing  poor  signal  quality  and  data  rate. 
Conversely,  because  of  the  crystalline  structure  of  the  1.5-ym  ^2^3 
particles ,  there  was  little  degradation  In  size  due  to  either  atomization  or 
thermal  effects. 

Once  It  was  determined  that  the  l.S-ym  particles  were  to  be  used,  a 
particle  neutralizer  was  employed  to  assess  the  effect  of  neutralizing  the 
charge  of  the  aerosol.  The  particle  neutralizer  In  the  present  tests  had  no 
significant  effect  on  the  size  of  the  generated  particle. 

Particle  Generator  Performance.  The  particle  generation  rate  was 
examined  by  utilizing  the  laser  system  shown  In  Figure  19b.  A  single  beam 
(green,  514.5  ym)  of  an  argon  Ion  laser  was  collimated,  focused,  and  directed 
through  the  combustor  test  section.  Light  scattered  by  particles  passing 
through  the  laser  beam  at  the  center  of  the  duct  was  focused  onto  a  0.25  mm- 
diameter  aperture  of  a  photomultiplier  tube  located  at  a  25°  angle  to  the 
laser  beam.  The  output  of  the  photomultiplier  was  amplified  and  fed  Into  a 
Schmitt  Trigger  which  provided  TTL  level  pulses  to  a  frequency  counter.  The 
number  of  particles  passing  though  the  area  of  sensltlvly  was  recorded  In  ten 
second  Intervals  and  scaled  by  the  area  ratio  (18,000:1)  of  the  test  section 
to  the  senslnt  area  to  obtain  a  total  particle  generation  rate. 

Particle  generation  rate  data  are  presented  In  Figure  19c  for  a 
period  of  time  exceeding  1  h.  Using  a  4  g/1  AI2O2-H2O  suspension,  25-pslg 
atomizer  pressure,  and  10:1  dilution  ratio,  the  mean  generation  rate  Is  1.1  x 
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10^  particles/ 8  with  rms  fluctuation  of  1.5Z. 

In  order  to  match  the  volumetric  concentration  of  Independently 
seeded  streams,  a  par£uneterlc  variation  on  the  atomization  seeding  system  was 
conducted.  The  particle  generation  rate  was  recorded  for  two  conditions:  (1) 
fixed  suspension  concentration  with  varying  differential  pressures  across  the 
atomizer  and  (2)  fixed  atomizer  differential  pressure  with  varying  suspension 
concentrations.  The  generation  rate  was  found  to  vary  linearly  with  both 
atomizer  differential  pressure  and  suspension  concentration  (Figure  19c)  for 
the  operating  conditions  considered. 

One  question  with  respect  to  this  technique  of  seed  generation  is 
the  amount  of  moisture  introduced  by  the  atomization  process.  At  a  10:1 
dilution  ratio,  the  relative  humidity  at  the  output  of  mixing  chamber  is  about 
50Z  at  72'’F  (0.8  weight  percent  of  water).  Since  the  aerosol  is  further 
diluted  by  an  additional  10:1  when  Introduced  into  each  stream  entering  the 
reactor,  the  presence  of  the  moisture  is  considered  negligible. 

Summary.  A  particle  generation  system  has  been  designed  and 
successfully  employed  to  Introduce  seed  in  uniform  concentration  to  flows 
typical  of  nonpremlxed  combustors  that  feature  multiple  stream  Injection.  An 
assessment  of  the  system  performance  resulted  in  the  following  observations: 

(1)  Particle  generation  rates  on  the  order  of  10^  particles/second  can 
be  produced  at  a  steady  rate  over  prolonged  periods  of  time. 

(2)  rms  fluctuatons  of  1.5Z  in  particle  generation  rate  are  typical  for 
operating  periods  in  excess  of  one  hour. 

(3)  Particle  generation  rate  is  linearly  dependent  on  (a)  atomization 
differential  pressure  for  fixed  suspension  concentration,  and  (b) 
suspension  concentration  for  a  fixed  atomizer  differential  pressure. 

(4)  Crystalline  rather  than  amophorus  aluminum  oxide  is  preferable  to 
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ensure  control  over  Che  final  size  of  Che  seed  particulate. 

Details  are  provided  In  Ikloka  (1982)  and  Ikloica,  Brum,  and  Samuelsen 
(1983). 

4.3.2  Dump  Combustor 

Introduction.  The  application  of  dump  combustor  configurations  In 
ramjets  has  prompted  Interest  In  the  fluid  dynamics  Involved.  As  part  of  this 
effort,  the  DCI  combustion  laboratory  has  performed  an  exploratory  analysis  of 
an  axlsymmetrlc  sudden  expansion  dump  combustor  utilizing  a  one-component 
laser  veloclmeter.  The  objectives  were  Co  (1)  develop  Insight  Into  Che 
aerodynamics  of  Che  flowfleld  for  both  non-reacting  and  reacting  flows,  and 
(2)  develop  a  database  for  the  purpose  of  model  evaluation  and  development. 
This  Cask  was  performed  under  a  special  funding  work  order  added  Co  the 
present  AFOSR  grant  by  the  Ramjet  Division  at  the  Air  Force  Wright 
Aeronautical  Laboratories  (AFWAL/PORT)  with  Dr.  F.D.  Stull  as  the  technical 
monitor. 

Experiment .  The  experimental  dump  combustor  Is  comprised  of  a  two- 
inch  stainless  steel  Inlet  section  which  dumps  Into  a  larger  quartz  test 
section  (see  Figure  20a) .  The  quartz  section  Is  removable  and  Is  capable  of 
accommodating  an  80  on  (3-lnch)  and  102  mm  (4-lnch)  diameter  tube.  The 
section  Is  91  cm  (3-feet)  long,  a  length  selected  to  contain  the  flame.  These 
systems  are  typically  operated  at  high  pressures  through  the  use  of  a  choked 
exit  nozzle.  Because  of  Che  air  flowrate  limitation  of  the  UCI  facility  and 
Che  safety  problems  Inherent  to  the  quartz  test  section,  Che  experiment  was 
performed  at  atmospheric  pressure,  thus,  negating  the  need  for  an  exit  nozzle. 

The  laser  veloclmetry  measurements  were  performed  utilizing  Che  UCI 
Combustion  Laboratory  laser  anemometer  system  (Figure  3b)  operating  In  the 
single  component  (axial)  mode. 
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The  experimental  test  configuration  and  teat  conditions  are 
tabulated  in  Table  V.  The  velocity  measurements  were  performed  at  17  axial 
stations  spaced  at  2  cm  Intervals  beginning  1  cm  downstream  of  the  dump 
plane.  At  each  axial  station,  radial  traverses  were  made  at  5  mm  Intervals. 
The  reacting  flow,  however,  was  limited  to  9  axial  stations  due  to  a  low  data 
rate. 

In  addition  to  the  low  data  rate  exhibited  by  the  hot  flow  run, 
radiative  heat  transfer  presented  experimental  challenges,  the  worst  being 
heat  transfer  to  the  photomultiplier  assembly.  This  was  alleviated  by 
increasing  the  focal  length  of  the  receiving  lows  to  250  na  from  120  mm, 
constructing  a  radiative  heat  shield  with  asbestos  backing  and  by  cooling  the 
optics  with  an  air  convection  system. 

Results .  Computer  generated  plots  of  axial  and  radial  velocity 
profiles  are  presented  in  Figure  20b  for  the  COLD  cases.  From  the  radial 
plots,  one  can  see  evidence  of  the  outer  recirculation  zones  associated  with 
the  sudden  expansion.  In  the  baseline  IS  m/s  case,  a  second  recirculation 
zone  is  apparent.  Also,  the  shear  layer  between  the  recirculation  zone  and 
the  main  flow  is  indicated  by  high  rms  velocity  levels.  Downstream  of  the 
recirculation,  the  velocity  profiles  approach  that  of  turbulent  plpeflow  with 
decreasing  values  of  RMS  velocity. 

The  recirculation  zone  in  the  core  of  HOT  flow  remains  unchanged 
(Figure  20c).  The  point  of  reattachment  occurs  at  approximately  7  step 
heights  downstream  of  the  dump  plane.  Note  the  acceleration  of  the  outer 
gases.  This  is  attributed  to  the  reactants  of  combustion  passing  normally 
through  the  conical  flame  front,  reacting,  and  expanding. 

Results  are  detailed  in  Roman  and  Samuelsen  (1982)  and  Roman  (1984). 


TABLE  V 


-  Dump  Combustor  Teat  Configurations  and  Test  Conditions  - 


CONFIGURATION 

<*1 

^2 

dj/dj 

Ai/A2 

L 

iiim(lnch) 

igm(lnch) 

mm(lnch) 

Baseline 

51(2) 

80(3) 

0.667 

0.44 

91(36) 

Off-Basellne 

51(2) 

102(4) 

0.500 

0.25 

91(36) 

MODE 

RUN 

VELOCITY 

m/s 

♦ 

CONFIGURATION 

COLD 

VI 

2.5 

NA 

Baseline 

V2 
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Baseline 

V3 

7.5 

NA 

Off-Basellne 

HOT 


HI 
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0.7 
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b)  Axial  and  Radial  Profiles:  COLD  Cases 
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Figure  20.  Dump  Conbustor  Results 
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c)  Axial  and  Radial  Profiles:  HOT  Case 
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Figure  20.  (concluded) 


4.3.3  NO  Sampling 

Introduction.  To  assess  the  conditions  under  which  transformations 
of  nitrogen  oxides  (NO  and  NO2)  may  occur  In  a  gaseous  sample  probe,  a 
controllable  laboratory  experiment  was  designed  and  used  In  an  earlier  study 
to  systematically  establish  the  conditions  under  which  NO  and  NO2  could  be 
transformed  and  the  extent  of  transformation  (Benson,  Samuelsen,  and  Peck, 
1976;  Samuelsen  and  Harman,  1977;  Benson  and  Samuelsen,  1977;  Samuelsen  and 
Benson,  1979).  Under  the  present  AFOSR  program  at  the  UCI  Combustion 
Laboratory  (Grant  78-3586),  the  Envlronics  Division  of  the  Air  Force 
Engineering  and  Services  Center  (AFESC/RDV)  provided  supplemental  funding  to 
modify  the  experiment  to  (1)  Include  water  vapor  In  the  simulated  combustion 
products,  (2)  provide  more  precise  control  over  the  temperature  and  history  of 
the  pollutant  mixture  prior  to  sampling,  and  (3)  simulate  gas  turbine  exhaust 
products  as  well  as  reacting  zone  gases,  automobile  exhausts,  and  boiler  flue 
gases. 

The  experiment  was  successfully  modified  and  operated.  Various 
compositions  of  simulated  exhaust  gases.  Including  those  representative  of  gas 
turbine  exhausts,  were  run  through  sample  probes  of  varying  materials  and 
analyzed  for  nitrogen  oxide  transformations.  This  analysis  was  performed  for 
various  probe  temperatures  over  a  range  (ISO'C  to  400'’C)  typical  of  the 
temperatures  experienced  by  probe  materials  In  practical  systems. 

Trials  were  run  using  2  m  (79-inche8)  long,  6.35  nm  (0.25-inch)  OD 
tubing,  and  performed  with  and  without  water  addition.  Quartz  and  316 
stainless  steel  were  the  principal  materials  evaluated  although  Inconel-600 
and  Nlckel-200  were  considered  as  well.  Type  316  stainless  steel  (316  SS)  Is 
the  material  commonly  used  for  sampling  purposes.  Tubing  stock  and  fittings 
are  readily  available,  316  SS  Is  durable  yet  easy  to  bend  during  fabrication, 


and  the  stainless  quality  Is  reslstent  to  corrosion.  However,  this  reslstence 
does  not  preclude  surface  reactions  that  could  alter  the  composition  of  a 
sample.  Quartz  tubing.  In  contrast.  Is  Inert  to  surface  attack,  and,  as  a 
result.  Is  attractive  for  sampling.  However,  the  vulnerability  of  quartz  to 
breakage  renders  the  material  Impractical  except  for  the  most  controlled 
conditions.  To  assess  whether  durable,  metallic  tubes  are  available  which  are 
less  reactive  than  316  SS  to  sample  composition  transformation,  two  nlckel- 
contalnlng  tube  materials  were  selected:  a  nickel  tube  (nlckel-200)  and  an 
alloy  of  nickel  (lnconel-600) .  An  Investigation  of  stainless  steel 
composition  (Samuelsen  and  Harman,  1977)  had  Indicated  that  the  nickel  content 
may  be  associated  with  an  Increase  In  resistivity  to  transformation. 

Experimental .  Carrier  gases  representative  of  boiler  and  automotive 
exhaust,  as  well  as  combustor  reaction  zones  are  herein  referred  to  as 
STANDARD  carriers;  those  representative  of  turbine  Jet  engines  are  referred  to 
as  the  carrier  TURBINE.  Compositions  and  concentrations  of  the  carriers, 
presented  In  Table  VI,  were  established  based  upon  surveys  of  practical 
devices  (Samuelsen  and  Harman,  1977;  Samuelsen,  Hamberg,  and  Osborn,  1978). 

All  STANDARD  carrier  gases  (exclusive  of  water  and  propylene)  were  mixed  and 
supplied  In  gas  cylinders  by  a  commercial  vendor.  Propylene  and  water.  If 
selected,  were  added  during  the  experiment.  The  TURBINE  carrier  gas  was  also 
supplied  by  a  commercial  vendor  except  for  the  water,  methanol,  and  propanol, 
each  of  which  was  added  In  the  conduct  of  the  experiment.  The  baseline 
composition  was  selected  based  on  exploratory  tests  on  the  condition  for  which 
NO  and  NO2  transformation  was  most  likely. 

The  experimental  system  Is  shown  In  Figure  21a.  The  experiment  Is 
performed  In  three  steps:  gas  preparation,  probe  exposure  (with  sampling 
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before  and  after  the  exposure),  and  analysis.  The  gas  preparation  begins  by 
selecting  the  appropriate  carrier  gas  by  means  of  a  toggle  switch  located  on 
the  control  panel.  Upstream  carrier  gas  pressures  are  maintained  constant  to 
generate  a  steady  flow  of  7  l/mln.  The  bulk  of  the  flow  Is  bubbled  through 
three  500  ml  fritted  glass  water  towers  to  ensure  saturation  If  water  addition 
("wet  trial”)  Is  to  be  performed.  The  temperature  of  the  water  In  the  towers 
Is  set  and  maintained  constant  per  each  trial  to  achieve  the  desired  water 
concentration.  For  a  dry  trial,  this  carrier  flow  Is  routed  around  the  water 
towers.  For  the  TURBINE  tests,  a  portion  of  the  carrier  Is  routed  through 
Erlenmeyer  flasks  (11)  used  to  house  liquid  methanol  and  propanol. 

Enroute  to  the  furnace  from  Che  water  Cowers ,  heating  Cape  Is 
wrapped  around  the  Teflon  tubing  to  prevent  premature  water  condensation.  In 
the  furnace  (Llndberg  Model  #51828),  the  carrier  gas  Is  circulated  until  It 
reaches  Che  set  furnace  temperature.  At  this  point,  Che  carrier  stream  enters 
a  mixing  venturi,  where  Che  pollutant  gases  (NO,  NO2  and  If  selected), 

are  Introduced  In  the  proper  amounts  from  hlgh-concentratlon  sources  by  means 
of  flow  restrictors.  A  sample  (31/ min)  Is  extracted  Immediately  after  mixing 
Is  completed.  (This  sample  location  Is  referred  to  as  SP2.)  A  flow 
of  41/  min  Is  run  chough  Che  probe  to  be  tested.  At  the  end  of  Che  probe, 
another  31/mln  sample  Is  extracted  at  sample  point  SP3,  leaving  ll/mln  to 
vent. 

Temperatures  are  monitored  at  SP2,  SP3,  In  Che  furnace,  and  Just 
downstream  of  Che  water  Cowers;  pressures  are  monitored  at  SP2,  SP3,  and  just 
upstream  of  the  water,  methanol /propanol  addition  system.  In  addition,  the 
carrier  gas  flow  rate  Is  monitored  Just  upstream  of  Che  water/methanol/propanol 


addition  system. 


TABLE  VI 


-  NO  Sampling  Carrier  Gas  Composition  - 


Carrier 

Composition 

CO2 

H,0 

HC 

NO 

NO2 

% 

% 

% 

ppm 

ppm 

ppm 

STANDARD 

12 

5 

0,  14 

0 

500 

75 

12 

5 

0,  14 

500  propylene 

500 

75* 

12 

0 

0,  14 

500  propylene 

500 

75 

12 

0 

0,  14 

0 

500 

75 

TURBINE 

4 

14 

0,  5 

18  methane 

0.5  butane 

34  ethylene 

12  propylene 
2.5  I-butene 

30  methanol 

500 

75 

3  propanol 

*  Baseline  composition  (0%  O2,  500  ppm  HC,  0%  water) 
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Figure  21.  NO  Sampling  Results 


b)  Sequential  Exposure  Cycles 
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Figure  21.  (concluded) 
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The  samples  from  SP2  and  SP3  are  immediately  cooled  by  water  Jackets 
In  a  specially  designed  water-cooled  quartz  extractor  to  60^C  and  transported 
through  the  furnace  wall  and  Into  a  Blue  H  oven  maintained  at  60 ”0.  The  oven 
houses  various  solenoids  as  well  as  a  carbon  NO2  — NO  converter.  Monitored 
temperatures  Include  that  of  the  oven,  the  cooling  water  (at  SP2,  SP3,  and  the 
mixing  venturi),  and  the  converter.  The  NO^  converter  jacket  temperature  Is 
controlled  by  a  constant  Input  voltage  from  a  varlac,  and  the  cooling  water 
flow  rate  Is  monitored  by  a  rotameter  just  upstream  of  the  Blue  M  oven. 

Instrumentation  used  to  analyze  the  samples  Includes  a  TECO  Model  10 
AR  Chemiluminescent  N0-N02-N0^  Analyzer,  a  Beckman  Model  951  NO/NO,^  Analyzer, 
and  a  Beckman  Model  590  HC/CO  Analyzer.  From  the  control  panel,  the  SP2  and 
SP3  samples  can  be  routed,  by  means  of  a  solenoid  off-on  switch,  to  the  TECO 
or  to  the  Beckman  951.  In  addition,  the  TECO  sample  can  be  routed  through  or 
around  the  converter.  By  the  use  of  metering  valves  and  rotameters  (located 
on  the  control  panel),  as  well  as  by  vents,  proper  flows  through  the  probe  and 
to  the  Instruments  are  maintained  regardless  of  the  sample  path.  For  the 
trials  In  which  propylene  Is  Introduced,  a  small  portion  of  the  Beckman  951 
sample  Is  extracted  and  routed  to  the  Beckman  590.  Water  Is  condensed  out  of 
the  samples  before  transportation  to  the  analyzing  Instruments  by  means  of  Ice 
baths . 

The  experiment  Is  run  In  the  following  manner:  with  the  appropriate 
carrier  gas  selected,  the  desired  concenteratlons,  If  any,  of  water  and 
methanol /propanol  are  added,  and  the  NO,  NO2  and  (If  selected)  are 

Introduced  at  55”C  to  the  carrier  stream  at  the  venturi  port.  The 
temperatures  of  the  cooling  water  (at  SF2,  SP3,  and  at  the  mixing  venturi), 
the  oven,  the  converter,  and  the  gas  just  downstream  of  the  water  towers  are 
noted.  The  furnace  Is  set  at  150*C  and,  when  this  temperature  Is  achieved, 
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the  temperature  of  the  furnace,  and  the  temperatures  and  pressures  of  SP2  and 
SP3  are  recorded.  The  following  data  are  collected:  a  TECO  analysis  of  the 
SP2  sample  with  the  converter  off ,  a  TECO  analysis  of  the  SP3  sample  with  the 
converter  off,  a  Beckman  951  analysis  of  the  SP2  sample,  a  TECO  analysis  of 
the  SP2  sample  with  the  converter  on,  and  a  TECO  analysis  of  the  SP3  sample 
with  the  converter  on.  The  furnace  temperature  Is  then  Incremented  50*C  (up 
to  400‘*C)  and  the  readings  above  are  repeated.  A.t  least  four  successive  probe 
trials  are  necessary  for  the  probe  response  to  stabilize. 

Results.  The  combination  of  sample  probe  materials  and  carrier  gas 
compositions  are  presented  In  Table  VII.  The  parameters  evaluated  were  probe 
exposure  cycles,  probe  preconditioning,  carrier  composition,  presence  of 
water,  and  probe  material. 

The  probes  tested  are  repeatedly  cycled  through  the  temperature 
range  described  above  until  a  steady-state  response  Is  achieved.  The  chemical 
transformation  observed  for  all  tests  conducted  Is  the  reduction  of  NO2  to 
NO.  "Reactivity"  Is  associated  with  the  change  of  NO/NO2  concentrations  from 
the  probe  Inlet  (SP2)  to  the  probe  outlet  (SP3).  If  NO2  reduction  Is 
experienced  In  a  particular  trial,  the  probe  Is  said  to  be  reactive,  with 
reactivity  varying  from  low  to  high  depending  on  the  extent  of  NO2  reduction 
and  the  temperature  at  which  It  first  occurs.  The  concentration  of  NO^  (NO  + 
NO2)  l9  also  monitored  and  found  to  be  conserved  In  all  the  tests  reported 
here;  although  In  previous  work,  as  noted  below,  destruction  of  NO^  was 
observed  at  the  elevated  test  temperatures. 

An  example  of  sequential  exposure  cycles  of  a  316  stainless  steel 
probe  Is  presented  In  Figure  21b  for  the  baseline  case  In  the  STANDARD 
matrix.  The  Initial  reactivity  of  the  probe  Is  very  high  (cycle  1)  with  an 
80Z  reduction  of  NO2  to  NO  and  the  first  evidence  of  change  at  250°C.  By  the 
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TABLE  VII  PROBE  TYPES  AND  CARRIER  COMPOSITIONS  TESTED 


Carrier 

STANDARD 

TURBINE 

Probe  Type 

Case 

Oj. 

Z 

HC,  ppm 

Water, 

Z 

Water,  Z 

0 

5 

0  500 

0 

14 

0  5 

316  SS,  or 

1 

• 

a 

a 

Ouartz 

2 

a 

a 

a 

3 

a 

a 

a 

4 

a 

a 

a 

5 

B* 

B 

B 

6 

a 

a 

a 

7 

a 

a 

a 

8 

a 

a 

a 

9 

a 

10 

a 

Inconel-600,  11 

B 

B 

B 

or  Nlckel-200 

a 


Baseline  composition 


fifth  cycle,  the  NO2  reduction  stabilizes  to  approximately  40Z,  and  the  first 
evidence  of  change  occurs  at  350°C.  This  behavior  is  typical  of  all  the  metal 
probes  tested.  Generally,  3  to  S  cycles  are  required  for  a  repeatable 
response  to  be  achieved.  In  contrast,  quartz  exhibits  no  change  from  that 
measured  for  the  first  cycle. 

For  the  results  reported  below,  the  probes  were  tested  as  received 
from  the  supplier.  However,  preliminary  tests  were  conducted  on  316  stainless 
steel  to  evaluate  the  effect  of  washing  the  probes  prior  to  testing.  The 
probes  were  washed  with  trlchlorylethylene,  rinsed  with  methanol,  and  dried  at 
room  temperature.  Probes  washed  In  this  fashion  proved  to  be  more  reactive. 
Extensive  discussions  with  manufacturers  of  stainless  steel  tubes  lead  to  the 
conclusion  that  this  washing  procedure  removes  residual  dirt  and  grease  from 
the  inner  tube  surface,  thereby  exposing  more  sites  to  contribute  to  the 
Initial  reactivity.  The  decision  to  use  probes  as  received  (rather  than 
washed)  for  the  principal  tests  was  predicated  on  the  likelihood  that  those 
probes  used  In  practical  systems  are  Installed  as  supplied. 

The  composition  of  the  carrier  gas  was  found  to  be  a  major  factor  In 
establishing  the  reactivity  of  the  316  stainless  steel.  The  results  for  the 
STANDARD  carrier  (Table  VI)  are  summarized  In  Figure  21c.  The  reactivity  Is 
highest  for  the  baseline  case,  lower  for  the  case  where  both  oxygen  and 
propylene  are  present,  and  Inconsequential  In  the  absence  of  the  propylene 
(for  both  0  or  5Z  oxygen).  In  contrast  to  the  316  stainless  steel,  the  quartz 
tube  Is  virtually  nonreactlve  for  all  STANDARD  carrier  compositions.  It  Is 
noteworthy  that  the  316  stainless  steel  and  quartz  reactivity  Is  less  than 
that  observed  and  reported  from  the  earlier  experiment  (Samuelsen  and  Benson, 
1979),  an  example  of  which  Is  presented  In  Figure  21d.  Two  factors  lead  to 
the  higher  reactivity  In  the  previous  work.  First,  the  temperature  was  not 


controlled  as  well  as  in  the  present  experiment.  The  probe  was  heated  In  a 
Llndberg  tube  furnace  rather  than  in  the  current  box  furnace.  As  a  result, 
the  probe  experienced  uneven  heating.  An  investigation  of  the  heating 
variation  revealed  that,  at  a  400 "C  set  point,  the  probe  could  locally  heat  to 
600“C  at  locations  where  the  tube  was  in  the  vicinity  of  the  heated 
filament.  The  design  of  the  present  experiment  ensured  that  the  temperature 
of  the  probe  was  maintained  Isothermal  at  the  set  point.  Second,  the  sample 
extracted  at  SP2  and  SP3  was  not  Immediately  cooled  In  the  previous  experiment 
and  reaction  could  persist  In  the  sampling  tube.  In  the  present  experiment, 
the  design  of  the  sample  ports  Included  water-cooling  to  Immediately  quench 
the  reaction. 

The  results  for  the  TURBINE  carrier  are  presented  In  Figure  21e. 

The  results  are  similar  to  the  STANDARD  carrier  where  both  oxygen  and 
hydrocarbon  are  present  albeit  In  differing  amounts.  The  first  transformation 
observed  occurs  at  300®C  and  reaches  a  45Z  reduction  of  NO2  at  400“C  in  316 
stainless  steel. 

The  results  reported  above  were  obtained  dry  (l.e.,  0%  H2O).  To 
assess  the  effect  of  water,  trials  were  also  conducted  with  the  water 
concentrations  listed  In  Table  VI,  14%  H2O  for  the  STANDARD  carriers,  and  5% 
H2O  for  the  TURBINE  carrier.  The  presence  of  water  had  no  significant  Impact 
on  reactivity  for  the  STANDARD  carriers  In  316  stainless  steel  except  for  the 
case  where  both  oxygen  (5%)  and  hydrocarbons  (500  ppm  C^Hg)  are  present. 

Here,  the  presence  of  water  substantially  Increased  the  reactivity  (Figure 
21f),  from  20%  to  85%  reduction  of  NO2  at  400**C.  The  Impact  of  water  In  the 


316  stainless  steel  on  the  TURBINE  carrier,  which  has  a  higher  oxygen  (14%) 
and  lower  hydrocarbon  (~  80  ppm),  was  not  significant.  In  quartz,  the  water 
had  no  Impact.  The  quartz  Is  virtually  passive  for  all  carriers  tested. 


Various  metallic  materials  were  tested  in  addition  to  the  316 
stainless  steel  to  determine  if  available  alternatives  to  quartz  would 
suppress  the  reduction  of  NO2  at  elevated  temperatures.  Both  nickel  and 
Inconel  tubes  were  tested  with  negative  results  (Figure  21g).  Both  materials 
exhibited  substantially  higher  reactivity  than  the  316  stainless  steel. 

Conclusions.  NO^  is  conserved  in  316  stainless  steel  and  quartz  tubes  up 
to  400®C.  Above  AQO^C,  can  be  destroyed. 

The  ratio  of  NO2/NO  la  conserved  up  to  300'’C  in  316  stainless  steel  and 
400'’C  in  quartz.  The  present  results  were  not  changed  appreciably  by  the 
presence  or  absence  of  water  with  one  notable  exception.  In  the  stainless 
steel,  the  reactivity  of  the  probe  (reduction  of  NO2  to  NO)  was  significantly 
Increased  above  300®C  (20-to  85-percent  at  400®C)  when  both  oxygen  and 
hydrocarbons  were  present  in  the  STANDARD  CARRIER. 

Stainless  steel  Is  more  rugged  than  quartz  and,  as  a  result,  more 
commonly  used  as  a  sampling  probe  material.  The  exploration  of  alternative 
metallic  tube  materials,  such  as  nickel  and  inconel,  was  not  successful  in 
identifying  a  durable  tube  material  that  was  less  reactive  than  316  SS.  Both 
nickel  and  Inconel  were  more  reactive. 

Based  on  the  results  of  this  study,  quartz  is  the  most  desirable  material 
to  use  for  sampling  at  elevated  temperatures.  Quartz  was  inert  to  all  the 
sample  compositions  and  sample  temperatures  explored  in  the  present 
experiment.  The  one  practical  exception  to  the  use  of  quartz  is  the 
vulnerability  to  breakage.  Stainless  steel  tubes  with  quartz  glass  linings, 
available  commercially,  represent  a  reasonable  solution  where  sample  integrity 
is  a  major  objective  and  reasonable  handling  of  the  probe  can  be  maintained. 

Should  quartz  or  quartz-lined  stainless  steel  not  be  an  option,  then  type 
316  stainless  steel  is  a  reasonable  alternative.  However,  if  stainless  steel 
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ts  used  at  temperatures  exceeding  200® C,  the  tube  should  be  water-cooled  to 
ensure  maintenance  of  sample  Integrity.  Care  must  be  excerclsed  to  cool  the 
tube  with  hot  water  (>  60®C)  to  ensure  that  water  In  the  combustion  products 
does  not  condense. 


SECTION  5 


tr 


SUMMARY 


The  development  of  combustors  that  are  both  fuel  efficient  and  fuel 
flexible  requires  spatially  resolved  measurements  of  velocity,  temperature, 
and  concentration  in  complex  flows  representative  of  both  premixed  and  non- 
premixed  systems.  Such  data  are  needed  to  provide  the  physical  Insight 
necessary  to  understand  the  basic  processes  of  turbulent  mixing  and  transport, 
and  support  the  evolution  of  modeling.  The  goal  of  the  AFOSR  program  at  the 
UCI  Combustion  Laboratory  was  to  develop  laboratory  model  combustors  and 
experimental  methodology  suitable  for  the  acquisition  of  the  desired 
information.  The  objectives  of  the  program  were: 

(1)  To  establish  and  evaluate  laboratory  model  combustors  for  both 
premixed  and  non-premlxed  fuel/air  injection  suitable  for  studies 
of  the  complex  flows, 

(2)  To  acquire  spatially  resolved  flowfield  data  using  state-of-the- 
art  diagnostics  for  the  purposes  of  model  verification  and 
providing  insight  into  the  transport  of  mass,  momentum,  and  energy 
in  complex  flows, 

(3)  To  conduct  supplementary  studies  to  support  the  development  of 
methodologies  required  for  measurements  in  complex  flows. 

Four  combustion  modules  were  considered:  the  Opposed  Jet  Combustor  (OJC) 
for  premixed  injection  and,  for  non-premixed  injection,  the  Centerbody 
Combustor  (CBC),  the  Centerbody  Combustor  with  swirl,  and  the  Dilute  Swirl 
Combustor  (DSC).  An  evaluation  of  combustor  performance  and  compatablllty 
with  a  list  of  criteria  for  laboratory  complex  flow  combustors  established  the 
DSC  as  a  viable  configuration  for  non-premlxed,  complex  flow  combustor  studies 
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in  which  swirl  was  a  major  aerodynamic  feature. 

Spatially  resolved  flowfleld  data  were  obtained  In  each  of  the  premlxed 

and  non-premlxed  modules.  Although  laser  anemometry  and  flow  vlsunllzatlon 

were  the  principal  diagnostic  tools  employed,  flow  mappings  of  species 

concentration  and  temperature  were  acquired  In  the  OJC  for  a  wide  range  of 

operating  conditions.  Studies  In  a  hydrod3rnaralc  (water-channel)  facility  were 

also  completed  for  both  the  OJC  and  CBC  geometries.  Extensive  two-component 

laser  anemometry  data  were  obtained  as  well  In  the  DSC  Including  measurement 

* 

of  the  Reynolds  stress.  In  a  complementary  program  with  the  National  Science 
Foundation  (NSF),  the  two-color  measurements  were  combined  with  a  fine-wire 
thermocouple  probe  to  yield  direct  measurement  of  both  axial  and  azimuthal 
heat  flux  In  the  DSC. 

Three  supplemental  studies  were  completed:  LA  Seeder,  Dump  Combustor, 
and  NO  Sampling.  The  LA  Seeder  study  addressed  the  development  and  evaluation 
of  a  seed  generator  for  laser  anemometer  measurements  In  a  reacting  flow  with 
two  or  more  separate  streams  of  fuel  and  air.  Such  flows  require  refractory 
particles  Injected  unlformally  In  each  of  the  Individual  streams  to  avoid 
biasing.  The  dump  combustor  study  was  an  Independent  element  to  the  AFOSR 
program  provided  by  the  Ramjet  Division  of  the  Air  Force  Wright  Aeronautical 
Laboratories  (AFWAL/PORT).  Single  component  laser  anemometry  data  were 
obtained  In  a  dump  combustor  configuration  with  two  area  ratios  for  both  non¬ 
reacting  and  reacting  flows.  The  NO  Sampling  study,  supported  by  the 
Envlronlcs  Division  of  the  Air  Force  Engineering  and  Services  Center 
(AFESC/ROV),  explored  the  chemical  transformations  of  nitrogen  oxides  that  can 
occur  while  sampling  combustion  products.  The  present  effort  expanded 
previous  work  conducted  for  AFOSR  to  (1)  conduct  experiments  with  more 
accurate  temperature  control,  (2)  consider  the  presence  of  water  vapor  In  the 
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simulated  combustion  products  as  well  as  dry  samples,  and  (3)  evaluate 
combustion  product  composition  representative  of  gas  turbine  exhausts. 
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